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ABSTRACT 


This  research  concerns  the  continuing  development  of  Seaweb  underwater  net¬ 
working.  In  this  type  of  wireless  network  the  radio  channel  is  replaced  by  an  underwater 
acoustic  channel  which  is  strongly  dependent  on  the  physical  properties  of  the  ocean  me¬ 
dium  and  its  boundaries,  the  link  geometry  and  the  ambient  noise.  Traditional  acoustic 
communications  have  involved  a  priori  matching  of  the  signaling  parameters  (e.g.,  fre¬ 
quency  band,  source  level,  modulation  type,  coding  pulse  length)  to  the  expected  charac¬ 
teristics  of  the  channel.  To  achieve  more  robust  communications  among  the  nodes  of  the 
acoustic  network,  as  well  as  high  quality  of  service,  it  is  necessary  to  develop  a  type  of 
adaptive  modulation  in  the  acoustic  network.  Part  of  this  process  involves  estimating  the 
channel  scattering  function  in  terms  of  impulse  response,  the  Doppler  effects,  and  the  link 
margin.  That  is  possible  with  the  use  of  a  known  probe  signal  for  analyzing  the  response 
of  the  channel.  The  estimated  channel  scattering  function  can  indicate  the  optimum  sig¬ 
naling  parameters  for  the  link  (adaptive  modulation).  This  approach  is  also  effective  for 
time  varying  channels,  including  links  between  mobile  nodes  (e.g.  two  submarines),  since 
the  channel  characteristics  can  be  updated  each  time  we  send  a  probe  signal. 
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EXECUTIVE  SUMMARY 


Seaweb  is  an  organized  network  of  battery-operated  acoustic  modem  nodes  de¬ 
ployed  on  the  seabed.  Those  nodes  support  bidirectional  communications  between  them, 
as  well  as  with  a  gateway  node.  Seaweb  is  designed  to  provide  command  and  control  to 
Unmanned  Underwater  Vehicles  (UUVs)  from  shore  facilities  or  surface  ships,  provide 
communications  between  submerged  submarines  and  land  bases,  and  enable  wide  area 
undersea  surveillance  in  littoral  waters. 

During  the  last  decade  numerous  experiments  took  place  in  many  different  acous¬ 
tic  channels.  One  very  interesting  result  from  those  experiments  is  that  communication 
performance  exhibits  time-varying  characteristics  strongly  dependent  on  the  variability  of 
the  underwater  channels  themselves.  In  this  thesis  we  analyze  the  characteristics  and  pa¬ 
rameters  of  the  underwater  channel  that  causes  the  communication  between  two  nodes  to 
be  so  environment  dependent. 

In  Chapter  II,  we  describe  how  the  sound  propagates  in  the  sea  and  consider  the 
colored  ambient  noise  existing  in  the  acoustic  medium.  We  analyze  how  a  signal  passing 
through  this  underwater  channel  is  distorted  both  in  the  time  and  frequency  domains.  We 
investigate  the  importance  of  the  impulse  response,  develop  the  characteristic  parameters 
of  the  channel  and  use  a  theoretical  model  to  describe  it. 

In  Chapter  III,  we  describe  the  Request  to  send  /  Clear  to  send  and  shake  protocol 
which  precedes  communication  in  Seaweb.  We  explain  how  we  can  incorporate  a  probe 
signal,  which  is  a  known  signal  of  special  format,  for  purposes  of  obtaining  the  channel 
parameters.  We  analyze  what  characteristics  a  signal  needs  in  order  to  be  used  as  a  probe 
signal,  and  then  refer  to  the  most  usual  ones.  We  describe  an  efficient  method  that  en¬ 
ables  estimation  of  the  scattering  function  of  the  channel.  This  method  uses  a  Direct  Se¬ 
quence  spread  spectrum  signal  as  a  probe  signal. 

In  Chapter  IV,  we  develop  an  artificial  channel  with  known  characteristics.  By 
passing  a  DSSS  signal  in  baseband  through  this  channel,  we  implement  the  previously 
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mentioned  method  to  the  received  signal  to  get  back  the  scattering  function  of  the  chan¬ 
nel.  The  scattering  function  explains  by  itself  every  aspect  of  the  channel.  We  also  exam¬ 
ine  how  the  addition  of  white  Gaussian  noise  influences  the  accuracy  of  the  estimation. 

In  Chapter  V,  we  analyze  data  from  a  Seaweb  experiment  (New  England  Shelf 
experiment  -  17-20  April  2000).  After  a  general  description  of  the  experiment  and  of  the 
various  probe  signals  sent,  we  estimate  the  characteristics  of  the  real  underwater  channel 
using  various  methods,  including  the  one  we  described  previously.  We  can  compare  then 
the  results  received  from  each  different  method.  The  data  were  recorded  in  ten  different 
ranges,  so  the  channel  estimation  takes  place  ten  times.  It  is  also  very  interesting  to  ob¬ 
serve  the  way  in  which  the  signal  is  affected  by  the  channel  in  gradually  increased  dis¬ 
tances.  Some  useful  conclusions  are  derived. 


I. 


INTRODUCTION 


A.  UNDERWATER  ACOUSTIC  NETWORKS 

Over  the  last  decade,  the  U.S.  Navy  has  begun  developing  underwater  acoustic 
networks  [1,2].  Those  networks  need  to  be  designed  to  provide  command  and  control  to 
Unmanned  Underwater  Vehicles  (UUVs)  from  shore  facilities  or  surface  ships,  support 
communications  between  submerged  submarines  and  land  bases,  and  enable  wide  area 
undersea  surveillance  in  littoral  waters.  The  product  of  this  development  is  Seaweb.  Sea- 
web  is  an  organized  network  of  battery-operated  acoustic  modem  nodes  deployed  on  the 
seabed  (Figure  1).  Those  nodes  support  bidirectional  communications  between  them,  as 
well  as  with  gateway  nodes.  The  nodes  are  networked  to  allow  the  hopping  of  data  from  a 
source  node  to  a  destination  node  through  a  combination  of  other  intermediate  nodes. 
During  the  last  decade  numerous  experiments  took  place  in  many  different  acoustic 
channels  [3].  One  very  interesting  aspect  of  those  experiments  is  that  communication  per¬ 
formance  exhibits  variability  related  to  the  variability  of  the  underwater  channels  them¬ 
selves.  On  those  experiments  a  network  of  nodes  cover  wide  areas  (such  as  a  five  by  fif¬ 
teen  nautical  miles  area)  exchanging  data  between  them. 


Figure  1.  Seaweb  illustration  (After  Ref.  2.) 
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B.  ADAPTIVE  MODULATION 

The  modem  presently  used  in  Seaweb  is  the  Benthos  modem  employing  a  modu¬ 
lation  scheme  involving  non-coherently  processed  frequency  shift  keying  (FSK)  [2,4], 
The  performance  of  any  other  modulation  scheme  has  to  be  compared  with  the  perform¬ 
ance  of  non-coherent  FSK,  since  it  provides  robust  communications  and  relatively  high 
data  rates.  The  performance  of  phase  coherent  signaling  degrades  more  rapidly  under  se¬ 
vere  channel  distortion  whereas  the  non-coherent  FSK  is  much  more  robust  [5,6],  In  the 
current  implementation  of  Seaweb,  the  parameters  of  the  communication  link  such  as 
frequency  band  of  operation,  modulation  scheme,  modem  output  power,  and  error  correc¬ 
tion  coding  type  are  determined  prior  to  deployment.  The  a  priori  choice  of  signaling  pa¬ 
rameters  tends  to  be  overly  conservative  and  non-optimal.  The  variability  of  the  acoustic 
channel  suggests  that  a  more  appropriate  way  to  deal  with  those  communication  parame¬ 
ters  is  to  determine  dynamically  which  combination  of  those  will  give  us  the  optimum 
communication  scheme.  This  means  that  the  signaling  parameters  can  change  for  each 
exchange  of  data  between  two  nodes,  depending  on  the  existing  characteristics  of  the 
channel.  As  a  result,  the  link  would  use  the  communication  scheme  with  the  highest  pos¬ 
sible  data  rate  and  the  minimum  possible  probability  of  error  (probably  on  the  order  of 
10“5).  In  the  same  view,  the  link  would  use  the  most  appropriate  frequency  band  and 
transmit  the  minimum  required  amount  of  power.  The  literature  refers  to  this  technique  as 
adaptive  modulation  [7,8].  The  proposed  scheme  of  adaptive  modulation  process  for 
Seaweb  is  presented  in  Figure  2.  The  parts  highlighted  with  yellow  are  analyzed  in  this 
thesis. 
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Figure  2.  Adaptive  modulation  Process 


C.  SCOPE  OF  THE  THESIS 

A  smart  modem  is  a  modem  that,  after  receiving  some  known  probe  signal  of  a 
special  format,  will  determine  the  characteristics  of  the  channel,  such  as  impulse  re¬ 
sponse,  Doppler  shift,  and  signal-to-noise  ratio.  Based  on  the  estimated  channel  condi¬ 
tions,  it  selects  the  communication  parameters.  The  goal  of  this  thesis  was  to  provide  an 
understanding  of  underwater  channel  estimation  for  determining  the  characteristics  if  the 
channel.  The  thesis  is  organized  as  follows. 

In  Chapter  II,  we  investigate  the  underwater  channel.  We  examine  the  factors  that 
make  this  channel  so  interesting  but  difficult  for  acoustic  communications.  We  compare 
the  underwater  channel  with  the  traditional  radio  channel  for  cellular  communications,  to 
improve  our  understanding. 

In  Chapter  III,  we  examine  the  method  of  channel  estimation.  We  describe  a  prac¬ 
tical  mechanism  for  obtaining  the  channel  characteristics  in  the  acoustic  modem.  We  con¬ 
sider  the  theory  of  channel  estimation. 

In  Chapter  IV  we  create  an  artificial  channel  with  known  characteristics  and  pass 
a  probe  signal  through  this  channel.  Then,  by  processing  the  received  signal,  we  obtain 
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characteristics  of  the  artificial  channel.  This  artificial  channel  is  a  test  channel,  useful  for 
confirming  our  implementation  of  the  channel  estimation  algorithm. 

In  Chapter  V,  we  process  data  from  an  actual  experiment  (New  England  Shelf  ex¬ 
periment  -  17-20  April  2000).  Using  the  code  developed  in  Chapter  IV,  we  determine  the 
characteristics  of  the  channel.  We  describe  the  conditions  of  this  experiment,  the  format 
of  the  data  used  and  also  refer  in  the  problems  encountered  in  extracting  the  scattering 
function  from  the  set  of  real  data. 

Chapter  VI  presents  a  summary  with  conclusions  and  goals  achieved.  Finally,  we 
discuss  the  goal  of  adaptive  modulation. 

In  the  Appendix  the  Matlab  codes  used  to  generate  the  plots  for  the  simulated  and 
the  real  channel,  are  presented. 
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II.  UNDERWATER  CHANNEL 


This  chapter  analyzes  the  basic  characteristics  of  the  underwater  channel.  Starting 
from  basic  acoustics  and  ray  propagation,  we  move  on  to  the  acoustic  communication 
channel  and  discuss  the  signal  distortion  effects  in  time  and  in  the  frequency  domain.  We 
also  define  the  various  types  of  fading  channels,  and  illustrate  the  most  usual  mathemati¬ 
cal  models  describing  their  behavior.  Then  we  define  the  impulse  response  and  demon¬ 
strate  how  it  influences  the  received  signal. 

A.  SOUND  PROPAGATION  IN  THE  OCEAN 

The  ocean  is  an  acoustic  waveguide  limited  above  by  the  sea  surface  and  below 
by  the  seafloor.  The  sea  surface  can  be  modeled  as  a  pressure  release  boundary  and  the 
sea  floor  as  a  second  fluid  medium  (Pekeris  Waveguide)  [9],  The  result  is  no  loss  of 
acoustic  energy  to  air,  but  almost  always  a  loss  of  energy  to  the  second  medium.  The  de¬ 
gree  of  this  effect  depends  on  the  characteristics  of  the  bottom  (sound  speed  -  density) 
and  on  the  incident  angle  of  the  acoustic  wave  to  the  bottom.  So  each  time  the  acoustic 
wave  impinges  the  bottom,  it  suffers  a  loss  in  strength. 

The  propagation  of  sound  in  the  ocean  can  be  described  in  various  ways,  but  for 
our  purposes  in  this  thesis,  we  follow  ray  theory,  an  approach  borrowed  from  optics.  This 
theory  is  based  on  the  assumption  that  energy  travels  along  reasonably  well  defined  paths 
through  the  medium.  However,  rays  are  not  exact  representations  of  waves,  but  only  ap¬ 
proximations  that  are  valid  under  certain  rather  restrictive  conditions  [10],  The  ocean  wa¬ 
ter  column  is  not  a  homogeneous  medium  and  the  sound  speed  varies  significantly  as 
function  of  depth.  So  the  acoustic  “rays”  in  the  ocean  refract  according  to  Snell’s  law. 
Snell’s  law  provides  a  simple  formula  for  calculating  the  ray  declination  angle  at  any 
depth  z  based  only  on  the  declination  angle  at  any  other  depth  and  knowledge  of  the 
sound  speed: 

c(zo)cos(0(z))  =  c(z)cos(0(zo))  (2.1) 

where  c(z )  is  the  sound  speed  at  depth  z.  A  general  rule  for  ray  propagation  is  that  a  ray 
always  bends  toward  the  neighboring  region  of  lower  sound  speed.  If  we  know  the  sound 
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speed  profile  of  the  water  column,  we  can  model  the  sound  propagation  in  this  environ¬ 
ment. 

An  example  of  ray  tracing  is  presented  in  Figure  3.  The  example  involves  sound 
propagation  in  a  120-meter  deep  sea.  The  transmitter  is  located  at  30  meters  depth,  and 
the  receiver  at  73  meters  depth  at  a  distance  12,800  meters  from  the  transmitter.  In  the 
left  plot,  the  sound  speed  profile  is  illustrated.  In  the  right  plot,  the  dominant  eigenrays 
are  shown.  The  eigenrays  are  the  rays  that  start  from  the  transmitter  and  pass  though  the 
position  of  the  receiver.  In  this  case  we  have  three  dominant  eigenrays.  Each  of  them  ex¬ 
hibits  different  values  of  attenuation,  phase  shift  and  propagation  delay.  The  red  eigenray 
arrives  at  the  receiver  first  and  experiences  the  least  attenuation.  Its  phase  shift  will  be 
close  to  zero.  Approximately  0.11  milliseconds  later  the  yellow  eigenray  reaches  the  re¬ 
ceiver,  with  a  phase  shift  close  to  pi  radians.  The  green  eigenray  reaches  the  receiver  0.55 
milliseconds  after  the  red  one,  with  a  phase  shift  close  to  n  as  well. 


Figure  3.  Ray  Tracing  example,  for  source  depth  at  30  meters 


The  intensity  of  the  propagating  sound  in  the  ocean  attenuates  as  a  function  of  the 
distance  r  between  the  transmitter  and  the  receiver.  The  term  that  describes  this  attenua- 
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tion  is  the  transmission  loss  (TL).  It  is  determined  by  combining  the  signal  loss  from 
source  to  receiver  due  to  a  combination  of  geometric  spreading  and  sound  absorption. 

Geometric  spreading  from  the  transmitter  is  in  the  form  of  spherical  spreading  up 
to  a  range  equal  to  the  water  depth  of  the  channel.  Spherical  spreading  loss  is  propor¬ 
tional  to  1 Jr1  .  Beyond  this  range,  cylindrical  spreading  approximates  the  propagation. 
Cylindrical  spreading  loss  is  proportional  to  1/r .  [11] 

The  absorption  of  sound  in  seawater  depends  on  numerous  parameters  such  as, 
temperature,  salinity,  depth,  pH,  frequency.  A  general  rule  is  that  as  frequency  increases 
the  absorption  of  sound  in  the  ocean  is  stronger.  Figure  4  illustrates  the  absorption  coeffi¬ 
cient  as  a  function  of  frequency  for  T  =  0°  C  and  T  =  20°  C,  for  a  depth  of  0  m,  pH  =  8 
and  S  =  35  ppt  (parts  per  thousand).  [9] 


Figure  4.  Absorption  in  seawater  -  Solid  line  is  for  T  =  0°  and  dashed  line  for  T  =  20° 

(After  Ref.  9.) 
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The  spectral  bandwidth  of  an  acoustic  communication  link  has  restrictions  on  its 
size.  Those  restrictions  have  two  origins.  With  higher  acoustic  frequencies  the  absorption 
increases  (see  Figure  4)  and  the  effective  distance  of  the  communication  link  will  be  very 
limited,  so  we  have  a  restriction  on  the  upper  side  of  the  band.  The  second  problem  is  the 
non-uniform  frequency  response  of  underwater  sound  projectors.  We  would  like  to  have 
a  relatively  flat  frequency  response  along  a  wide  frequency  band.  This  is  usually  hopeless 
because  of  the  way  those  transducers  are  built.  At  typical  communications  frequencies 
between  15  and  30  kHz,  they  can  have  a  relatively  flat  response  along  a  bandwidth  of  15 
kHz,  at  best.  An  example  of  a  typical  frequency  non-uniform  response  of  a  Tonpilz  trans¬ 
ducer  with  a  matching  layer  is  shown  in  Figure  5  [12]. 


p  (Hz) 

Figure  5.  Tonpilz  transducer  frequency  response  (After  Ref.  12.) 


B.  NOISE 

The  ambient  noise  of  the  ocean  is  not  Gaussian  and  colored.  Below  500  Hz  the 
major  contribution  to  ambient  noise  is  from  distant  shipping  and  biological  noise.  From 
500  Hz  to  50  kHz  the  local  sea  surface  is  the  strongest  source.  This  is  actually  the  fre- 
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quency  band  of  interest,  since  most  acoustic  modems  operate  in  this  band.  As  a  result,  the 
channel  noise  varies  depending  on  the  conditions  (sea  state,  wind),  as  illustrated  in  Figure 
6.  Above  50  kHz  the  ocean  turbulence  and  the  thermal  agitation  of  the  water  molecules 
are  the  predominant  noise  source.  [13] 


Figure  6.  Deep  water  ambient  noise  (After  Ref.  9.) 

Traditionally  the  theory  of  communications  assumes  Additive  White  Gaussian 
Noise  (AWGN).  However  AWGN  is  not  representative  of  the  acoustic  channel.  The 
noise  is  neither  white  nor  Gaussian.  The  result  is  that  communications  performance  is 
worst  and  also  the  analysis  is  much  more  complicated. 

C.  SIGNAL  DISTORTION  DUE  TO  MULTIPATH  PROPAGATION 

In  the  case  of  underwater  acoustic  communications,  the  signal  is  carried  by  acous¬ 
tic  pressure  waves.  In  this  channel  there  are  some  interesting  effects  taking  place,  which 
put  severe  limitations  in  the  communications.  The  first  effect  is  illustrated  in  Figure  3. 
We  see  that,  for  the  geometry  shown,  the  acoustic  waves  reach  the  receiver  following 
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three  different  paths.  Destructive  interference  by  the  multipath  propagation  structure  cre¬ 
ates  severe  fading  in  the  channel  [14,  15]. 

Fading  is  caused  by  interference  of  two  or  more  replicas  of  the  transmitted  signal 
arriving  at  the  receiver  at  slightly  different  times  following  different  paths.  These  differ¬ 
ent  components  are  called  multipaths.  They  combine  at  the  receiver  to  give  a  resultant 
signal,  which  can  vary  widely  in  amplitude  and  phase,  depending  on  the  distribution  of 
the  intensity  and  relative  propagation  time  of  the  waves  and  the  bandwidth  of  the  trans¬ 
mitted  signal.  The  multipath  signal  viewed  in  the  frequency  domain,  exhibites  different 
spectral  components  of  the  signal  being  affected  differently  by  the  channel.  In  other 
words  the  frequency  response  of  the  channel  is  not  flat  over  the  bandwidth  of  the  signal. 

Another  important  characteristic  of  the  multipath  propagation  is  the  time  variation 
in  the  structure  of  the  acoustic  medium  (waves,  wind,  current)  and  the  motion  of  either 
the  receiver  or  the  transmitter.  As  a  result,  the  signal  passing  through  the  underwater 
channel  is  distorted  both  in  the  time  and  frequency  domains. 

1.  Energy  Time  Spread 

Consider  the  simplified  case  of  Figure  3.  We  assume  a  sinusoidal  signal  with  a 
transmit  duration  of  one  millisecond.  The  acoustic  wave  (following  Figure  3)  follows 
three  different  paths.  Each  path  has  some  attenuation,  some  phase  shift  and  some  propa¬ 
gation  delay.  The  transmitted  pulse  is  illustrated  on  the  left  side  of  Figure  7  and  the  re¬ 
ceived  pulse  on  the  right  side.  As  we  can  see,  the  received  pulse  is  distorted  and  its  en¬ 
ergy  is  spread  in  time. 


10 


Figure  7.  Multipath  effect  on  a  sinusoidal  pulse 


2.  Doppler  Shift  -  Doppler  Spread 

In  addition  to  the  energy  time  spread,  other  important  phenomena  are  observed  in 
the  underwater  channel.  They  have  to  do  with  the  distortion  of  the  signal  in  the  frequency 
domain.  Their  origin  is  the  time  variations  in  the  structure  of  the  acoustic  medium  and  the 
motion  of  either  the  receiver  or  the  transmitter. 

The  motion  of  the  receiver  or  transmitter  gives  us  the  well-known  Doppler  shift 
effect  wherein  the  center  frequency  of  the  carrier  is  shifted.  The  shift  is  positive  if  re¬ 
ceiver  and  transmitter  are  coming  closer  to  each  other,  and  negative  if  they  are  moving 
away  from  each  other.  The  amount  of  shift  is  given  by  the  following  equation: 

fd  [Hz]  =  ^vcos(6>).  (2.2) 

It  depends  on  the  frequency/  the  speed  of  sound  c,  the  relative  speed  v  between  the  re¬ 
ceiver  and  transmitter  v ,  and  the  spatial  angle  6  between  the  direction  of  motion  and  the 
direction  of  arrival.  Since  for  each  path  the  wave  arrives  from  a  different  angle  (in  gen¬ 
eral),  there  is  a  different  Doppler  shift  associated  with  each  path. 

The  time  variations  of  the  structure  of  the  acoustic  medium  give  us  a  relative  phe¬ 
nomenon,  which  is  called  Doppler  spread.  Assume  that  we  send  a  pure  tone  through  the 
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channel.  If  the  channel  is  time  invariant,  we  do  not  notice  any  spectral  broadening  in  the 
received  tone.  However,  time  variations  of  a  channel  result  in  a  broadening  of  the  spectral 
line.  This  effect  is  what  we  call  Doppler  spread.  The  Doppler  spread,  just  like  the  Dop¬ 
pler  shift,  has  a  different  value  for  each  path.  [16] 

D.  IMPULSE  RESPONSE  PROFILE  -  IMPORTANCE 

As  we  discussed  earlier,  each  path  has  a  different  attenuation,  phase  shift  and  de¬ 
lay.  Therefore,  if  a  real  bandpass  signal  ,v(?)  is  sent  though  the  channel,  the  received  sig¬ 
nal  can  be  expressed  as 

rM  =  2X(*M  t~Tn(t)\  (2-3) 

n= 0 

In  this  equation,  N  is  the  total  number  of  paths  contained  in  the  underwater  channel, 
an  (?  )  is  the  attenuation  of  the  »-th  path,  which  is  a  function  of  time,  and  rn  (?)  is  the 

time  delay  associated  with  the  n- th  signal  path  and  is  a  function  of  time  as  well,  with 
r0  =  0  .  If  we  express  this  signal  in  the  baseband,  then  the  baseband  equivalent  received 
signal  has  the  form: 

N 

r{t)  =  Y}an{t)^(t-Tn{t))^\}{0)c^Tn{t)  +  (Pn{t))\  (2.4) 

n= 0 

where  s(t-rn  (?))  is  the  baseband  form  of  the  transmitted  signal  and  <pn  (?)  is  the  phase 

shift  of  the  n- th  path  and  is  a  function  of  time.  We  define  the  complex  baseband  impulse 
response  of  the  multipath  channel  as:  [4] 

Mr’0  =  El  an  (0 1  exp[/(^.rH  (t)  +  <pH  (?))]<? (r-r„  (?)).  (2.5) 

n= 0 

From  the  last  relation  we  can  see  that  the  impulse  response  of  the  channel  is  the  superpo¬ 
sition  of  the  impulse  responses  of  the  individual  paths.  This  is  a  very  important  quantity 
for  our  analysis  since,  after  determining  its  form,  we  are  able  to  derive  all  other  character¬ 
istic  parameters  of  the  channel  using  this  result.  More  specifically,  we  may  determine  the 
time  and  frequency  effects  of  the  underwater  multipath  fading  channel  on  our  signal. 
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E.  UNDERWATER  CHANNEL  CHARACTERISTICS  PARAMETERS 


The  underwater  channel  can  be  characterized  by  at  least  two  parameters.  The  first 
characterizes  the  time  variations  and  the  other  the  frequency  variations  of  the  channel. 

The  Fourier  transform  of  the  impulse  response  /i(r,t)  with  respect  to  the  time  de¬ 
lay,  is  given  by  the  relation 

+00 

H(f,t)  =  J  h(r ,t)e~l2icft dt.  (2.6) 

—oo 

The  autocorrelation  of  the  channel  impulse  response  Fourier  transform  (with  respect  to 
time  delay),  is  given  by  the  relation 

(2.7) 

where  E{Q  stand  for  the  expected  value  of  the  function  inside  the  brackets. 

If  we  take  the  inverse  Fourier  transform  of  SH(Af;  t,  tx )  with  respect  to  A/ ,  we  obtain 

Ph  (r;  f, 0  =  SH  (A/;  t, tx ) en*Afd  (A/).  (2.8) 

When  the  channel  is  wide  sense  stationary  in  the  t  variable,  Ph  (r;t,tj)  =  Ph  (r,At)  where 

At  =  t  - 1, .  In  the  case  when  At  =  0 ,  Ph  (r)  is  the  average  power  out  of  the  channel  as  a 

function  of  the  time  delay  r  .  It  is  called  the  power  delay  profile  or  multipath  intensity 

profile  [Ref.  4,  16].  The  range  of  r  over  which  the  power  delay  profile  is  nonzero  is  the 
multipath  spread  Tm  of  the  channel.  The  multipath  spread  is  the  first  of  the  two  important 

parameters  and  describes  the  time  dispersive  nature  of  the  channel.  Some  typical  values 
of  the  multipath  spread  for  an  underwater  channel  are  10-20  ms,  compared  to  the  cellular 
wireless  communication  channel,  where  they  are  1-10  ps  [3,  4],  The  greater  the  multipath 
spread,  the  more  dispersive  the  channel  will  be. 

Now  consider  SH  (A/;  At)  .  In  the  case  when  At  =  0  ,  SH  (A/ )  is  the  frequency 
correlation  function.  The  range  of  A f  over  which  the  frequency  correlation  function  is 
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greater  than  some  defined  value  is  the  coherence  bandwidth  of  the  channel  Bc .  Since 
SH  (A/ )  and  Ph  (t)  are  a  Fourier  transform  pair,  the  relation  Bc  « \/Tm  will  hold.  Conse¬ 
quently,  both  parameters  describe  the  time  dispersive  nature  of  the  channel.  The  coher¬ 
ence  bandwidth  is  a  measure  of  the  range  of  frequencies  over  which  the  channel  can  be 
considered  flat,  which  means  that  the  channel  will  pass  all  spectral  components  with  ap¬ 
proximately  equal  gain  and  linear  phase. 

Again  consider  SH  (A/;  At) .  In  the  case  when  A/  =  0 ,  SH  (At)  is  the  space-time 
correlation  function,  which  represents  the  correlation  of  a  single  sine  wave  with  itself 
over  time  At .  The  time  over  which  SH  (At)  is  essentially  unity  is  the  coherence  time  of 

the  channel  (At)  .  Coherence  time  is  the  parameter  which  describes  the  time- varying  na¬ 
ture  of  the  channel.  It  is  a  measure  of  the  time  duration  over  which  the  channel  attenua¬ 
tion  and  delay  are  essentially  constant,  so  the  received  amplitude  and  phase  are  constant 
over  a  period  of  (At)c  seconds.  If  we  take  the  Fourier  transform  of  SH  (A/;  At)  with  re¬ 
spect  to  At  we  obtain 

Ph  (A/;v)  =  [jlf  (A/;  At) e~i2m/Atd  (At ).  (2.9) 

In  the  case  when  A/  =  0 ,  Ph  (  v  )  is  the  Doppler  power  spectrum  of  the  channel.  The  range 
of  v  over  which  Ph  (v)  is  essentially  nonzero  is  the  Doppler  spread  Bd  of  the  channel. 
The  Doppler  spread,  as  we  discussed  before,  is  a  measure  of  the  spectral  broadening 
caused  by  the  time  rate  of  change  of  the  underwater  channel.  Since  SH  (At)  and  Ph  (v) 

are  Fourier  transform  pair,  the  relation  Bd  « l/(At)c  will  be  true.  As  a  result  of  that,  both 

parameters  describe  the  time-varying  nature  of  the  channel.  [4] 

F.  DIFFERENT  TYPES  OF  FADING  CHANNELS 

We  can  characterize  a  multipath  fading  channel.  This  can  be  done  by  comparison 
of  the  parameters  of  the  channel  to  the  characteristics  of  the  communication  signal,  spe¬ 
cifically,  the  signal’s  bandwidth  and  symbol  period. 
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Time  dispersion  due  to  multipath  causes  the  transmitted  signal  to  experience  ei¬ 
ther  flat  or  frequency  selective  fading.  If  the  channel  has  a  constant  gain  and  linear  phase 
response  over  a  bandwidth  which  is  greater  than  the  bandwidth  of  the  transmitted  signal 
(i.e.,  Bc  is  large  in  comparison  with  the  BW),  the  channel  is  said  to  be  ffequency- 
nonselective  or  flat  fading.  In  the  time  domain  this  means  that  all  of  the  multipath  com¬ 
ponents  arrive  within  the  symbol  duration.  On  the  other  hand,  if  B  is  small  in  compari¬ 
son  with  the  BW,  significant  distortion  of  the  signal  occurs  and  the  channel  is  said  to  be 
frequency-selective.  In  this  case  successive  pulses  interfere  with  each  other. 

Figure  8  illustrates  the  concept  of  flat  and  frequency  selective  fading. 


(ii)  Typical  frequency  selective  fading  case  (Bc  <  BW) 

Figure  8.  Examples  of  frequency  non-selective  /  selective  fading  channels 

Depending  on  how  rapidly  the  transmitted  signal  changes  as  compared  to  the  rate 
of  change  of  the  channel,  a  channel  may  be  classified  either  as  a  fast  fading  or  slow  fad¬ 
ing  channel.  If  the  symbol  duration  is  smaller  than  the  coherence  time,  then  the  received 
amplitude  and  phase  are  effectively  constant  for  the  duration  of  at  least  one  symbol  and 
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the  channel  is  said  to  be  slowly  fading.  But  if  the  received  amplitude  and  phase  fluctuate 
over  time  periods  that  are  short  compared  to  the  duration  of  a  symbol,  the  channel  is  said 
to  be  fast  fading. 

The  underwater  channel  is  always  strongly  frequency  selective,  with  the  large 
multipath  spread  likely  to  cause  Inter  Symbol  Interference  (depending  on  the  signaling 
rate)  up  to  several  tens  of  symbol  intervals.  The  channel  depending  on  the  conditions  and 
its  geometry  can  be  slow  or  fast  fading  with  the  worst-case  scenario  being  in  shallow  wa¬ 
ters  and  under  rough  weather.  [13] 

G.  POSSIBLE  MODEL  OF  UNDERWATER  CHANNEL 

Since  the  time  variations  of  the  multipath  channel  appear  to  be  unpredictable  to 
the  user,  it  is  reasonable  to  characterize  this  time  variant  channel  statistically.  When  the 
impulse  response  /i(r,t)  is  modeled  as  a  zero  mean  complex  valued  Gaussian  process 

then  the  envelope  of  |/i(r,t)|  at  any  instant  t  is  Rayleigh-distributed.  So  the  Rayleigh  dis¬ 
tributed  case  describes  a  channel  in  which  there  is  no  line-of-site  direct  path.  This  type  of 
channel  is  referred  as  a  Rayleigh  fading  channel.  If  there  are  fixed  scatterers  or  signal  re¬ 
flectors  in  the  medium  in  addition  to  randomly  moving  scatterers,  h(r,t )  can  no  longer 

be  modeled  as  having  zero  mean.  In  this  case  the  envelope  of  h(r,t )  at  any  instant  t  is 

Rayleigh-distributed,  and  the  channel  is  characterized  as  a  Ricean  fading  channel.  In  the 
Ricean  fading  channel,  we  have  better  communication  performance  than  in  the  Rayleigh 
case. 

At  first  glance,  the  underwater  channel  would  seem  to  be  Ricean  fading,  since  the 
sea  surface  and  seafloor  represent  fixed  scatterers  in  the  acoustic  medium.  However,  ex¬ 
periments  show  that  the  behavior  of  the  channel  depending  on  the  conditions  resembles 
either  a  Ricean  or  a  Rayleigh  channel.  In  deep  waters  and  when  the  wind  is  feeble,  the 
underwater  channel  behaves  as  a  Ricean  fading  channel.  On  the  other  hand,  in  shallow 
waters  with  strong  winds,  the  channel  behaves  as  a  Rayleigh  fading  channel.  The  second 
case  is  worst  and  the  performance  is  poorer.  [4] 
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H.  CHAPTER  SUMMARY 

In  this  chapter,  we  analyzed  the  underwater  channel  and  all  its  characteristics  that 
make  it  so  interesting  and  difficult  for  acoustic  communications.  Next,  we  examine  the 
method  of  channel  estimation  by  describing  a  practical  mechanism  for  obtaining  the 
channel  characteristics  in  the  acoustic  modem. 
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III.  CONTEXT  OF  CHANNEL  ESTIMATION 


In  this  chapter,  we  examine  the  method  of  channel  estimation.  We  describe  a 
practical  algorithm  for  obtaining  the  channel  characteristics  in  the  acoustic  modem. 

A.  RTS  /  CTS  PROCEDURE 

The  Seaweb  underwater  network  provides  for  an  RTS/CTS  handshaking  proce¬ 
dure  for  setting  up  data  transmissions  as  depicted  in  Figure  2.  This  procedure,  in  the  cur¬ 
rent  form  of  Seaweb,  is  implemented  in  as  an  exchange  of  9-byte  utility  packets.  Prior  to 
communication,  two  acoustic  modems  (nodes  in  the  Seaweb  network)  can  perform  this 
handshake  in  order  to  establish  a  link.  [17] 

Consider  a  node  A  which  intends  to  send  data  to  node  B.  First,  node  A  transmits  a 
Request  to  Send  (RTS)  utility  packet.  In  the  current  structure  of  the  acoustic  modem,  the 
RTS  message  is  converted  from  binary  to  M-ary  data  symbols,  and  is  passed  through  a 
convolutional  encoder  with  code  rate  Vi  and  an  interleaver,  which  scrambles  the  coded 
symbols  in  order  to  make  the  link  more  robust  against  fading  effects.  The  modulation 
used  is  MFSK.  A  synchronization/acquisition  signal  is  appended  at  the  beginning  of  the 
signal.  The  RTS  packet  passes  through  the  channel  and  is  distorted  due  to  multipath 
spread.  Since  this  signal  is  very  important  for  the  link,  we  have  to  be  sure  that  it  will 
reach  the  node  B  correctly.  For  that  reason,  we  use  a  long  symbol  period  (on  the  order  of 
50  milliseconds),  in  order  to  ensure  that  the  symbol  period  is  longer  than  the  multipath 
spread  of  the  underwater  channel.  Upon  reception,  the  signal  is  demodulated  by  use  of 
noncoherent  means,  and  is  deinterleaved  and  decoded.  After  reception,  node  B  acknowl¬ 
edges  back  that  it  is  ready  to  receive,  by  returning  a  Clear  to  Send  (CTS)  utility  packet. 
[18] 

This  process  is  what  is  called  RTS/CTS  handshaking.  The  RTS  message  is  used  to 
wake  up  node  B  and  prepare  it  for  the  reception  of  the  data.  The  9  bytes  conveys  informa¬ 
tion  about  the  data  packet,  and  other  housekeeping  data.  The  CTS  likewise  conveys  over¬ 
head  information.  Node  B  returns  the  chosen  communication  parameters  to  node  A  as  a 
specification  embedded  in  the  CTS  utility  packet.  For  the  purposes  of  this  thesis,  a  part  of 
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the  RTS  signal  is  used  for  channel  estimation  at  node  B.  In  the  future,  the  channel  esti¬ 
mate  is  the  input  information  to  the  acoustic  modem  for  determining  the  optimal  commu¬ 
nication  parameters. 

B.  PROBE  SIGNALS 

In  order  to  determine  the  channel  characteristics,  node  A  must  send  node  B  a  sig¬ 
nal  known  to  B.  This  signal  is  referred  to  as  a  probe  signal  and  must  have  a  special  for¬ 
mat.  The  probe  signal  passess  through  the  underwater  channel  and  it  is  distorted  in  time 
and  frequency  as  seen  in  Chapter  II.  Since  node  B  knows  in  advance  what  the  form  of  the 
probe  is,  it  can  determine  what  the  effect  of  the  channel  was  by  the  use  of  appropriate 
signal  processing. 

As  discussed  earlier,  the  RTS  signal  is  very  important  in  Seaweb,  so  the  informa¬ 
tion  carried  during  the  RTS  procedure  has  to  reach  the  target  node  without  errors.  Node  B 
demodulates  the  RTS  signal,  and  reconstructs  a  clean  replica  of  the  waveform  transmitted 
by  node  A.  This  waveform,  or  an  appended  special-purpose  waveform,  serves  as  the 
channel  probe. 

The  probe  signal  must  be  suitable  to  estimate  the  dynamics  of  the  underwater 
channel.  Wideband  probe  signals  provide  high  resolution  in  time  and  frequency  and  they 
are  often  used  in  practical  systems  to  measure  the  channel  characteristics.  Typical  wide¬ 
band  signals  for  this  purpose  can  be  a  Linear  Frequency  Modulated  (LFM)  chirp  or  a 
pseudorandom-noise-spread  signal  (such  as  a  Direct  Sequence  Spread  Spectrum  signal  - 
DSSS). 

1.  LFM  Chirp 

The  first  type  of  wideband  signal  that  can  be  used  as  a  probe  to  the  channel  is  the 
Linear  Frequency  Modulated  (LFM)  chirp,  a  sinusoidal  signal  with  frequency  sweeping 
with  time  in  a  linear  way.  The  LFM  chirp  signal  has  quadratic  phase.  The  form  of  an 
LFM  chirp  is 

x(t)  =  Acos^at2  +  f3t  +  y^.  (3.1) 

Its  instantaneous  frequency  is  f(t)  =  2at  +  j3.  We  can  see  that  the  frequency  changes 

linearly  with  time.  As  a  result  for  our  case,  the  overall  signal  is  wideband.  Both  the  fre- 
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quency  band  and  resolution  depend  on  the  values  of  the  time  duration  of  the  LFM  chirp 
and  on  the  parameter  a.  In  the  example  of  LFM  chirp  shown  in  Figure  9,  the  parameters 
are  set  so  that  the  sweeping  frequency  is  in  the  range  100  to  400  Hz  and  the  chirp  dura¬ 
tion  is  one  second.  [19] 


LFM  CHIRP 


FREQUENCY 


Figure  9.  Example  of  an  LFM  chirp 


2.  DSSS  Signal 

The  second  type  of  signal  that  can  be  used  for  channel  sounding  is  the  direct  se¬ 
quence  spread  spectrum  (DSSS)  signal.  This  type  of  signal  is  obtained  by  mixing  a  car¬ 
rier  signal  with  a  pseudonoise  (PN)  random  sequence. 

The  characteristics  of  the  PN  sequence  have  to  be  like  those  of  a  true  random  bi¬ 
nary  sequence.  The  most  important  figure  of  merit  for  the  PN  sequence  is  the  autocorrela¬ 
tion  function.  It  has  to  resemble  that  of  a  true  random  binary  sequence.  The  comparison 
of  the  two  autocorrelation  functions,  as  presented  in  Figure  10,  shows  great  similarity, 
except  for  the  periodicity.  Also,  the  greater  the  number  of  chips  N  in  the  PN  sequences, 
the  better.  A  PN  sequence  can  be  generated  by  the  use  of  an  n-stage  shift  register  where 
the  output  of  each  stage  are  properly  connected  or  not  connected  to  an  exclusive-or  gate 
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whose  output  is  fed  back  to  the  input  of  the  shift  register.  The  resulting  sequence  is  called 
a  maximal-length  sequence  or  just  a  m-sequence.  In  a  set  of  length-A  m-sequences,  some 
will  have  better  crosscorrelation  properties  and  those  are  called  preferred  m-sequences. 
Combining  appropriate  sets  of  those,  we  can  get  another  set  of  PN  sequences  which  are 
called  Gold  sequences.  [Ref.  4] 


(i)  Auto  correlation  function  of  a  polar  random  binary  wave 


(i)  Autocorrelation  function  for  a  PN  sequence  with  chip  duration  Tc 


Figure  10.  Comparison  of  the  autocorrelations  of  PN  sequence 
and  random  binary  sequence 

The  generation  of  a  DSSS  bandpass  signal  (in  the  time  domain)  is  now  briefly  ex¬ 
plained  following  Figure  11.  Assume  that  we  want  to  transmit  three  da  1,  0  and  1.  We 
transform  the  binary  data  in  the  polar  binary  wave  (1  and  -1)  as  shown  in  the  first  quar¬ 
ter  of  the  figure.  The  duration  of  the  bit  is  0.75  milliseconds,  which  corresponds  to  a 
bandpass  null-to-null  bandwidth  of  2.67  kHz.  In  the  second  part  of  the  figure  is  an  m- 
sequence  with  a  length  of  15  chips.  It  is  important  to  notice  that  15  chips  occupy  a  time 
duration  of  0.75  milliseconds  (i.e.,  the  duration  of  one  data  bit).  In  order  to  create  the 
baseband  DSSS  signal  we  multiply  the  two  binary  waveforms  shown  in  the  first  half  of 

the  figure;  the  resulting  signal  is  shown  on  the  third  line.  After  mixing  the  baseband 
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waveform  with  a  carrier  of  160  kHz,  we  get  the  DSSS/BPSK  form  which  is  shown  in  the 
last  quarter  of  the  figure.  Since  the  chip  duration  is  0.05  milliseconds,  which  is  15  times 
smaller  than  the  bit  duration,  the  null-to-null  bandwidth  of  the  resulting  DSSS/BPSK  sig¬ 
nal  is  1 5  times  larger  than  that  of  an  unspread  BPSK  signal,  equal  to  20  kHz.  In  this  ex¬ 
ample  the  length  of  the  spreading  code  is  relatively  small,  since  in  practice  we  use  much 
longer  PN  sequences.  For  acoustic  communications,  Gold  codes  of  length  2047  are  some¬ 
times  used. 


DATA  BUS 


0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8  2  2.2 

TIME  IN  msecs 

Figure  11.  Direct  Sequence  Spread  Spectrum  Signal  Generation 


In  order  to  demodulate  the  DSSS/BPSK  waveform  it  first  has  to  be  despread.  The 
receiver’s  knowledge  about  the  chipping  signal  is  necessary.  If  we  know  the  chipping  se 
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quence,  the  spread  signal  is  multiplied  with  the  aligned  and  synchronized  PN  sequence 
and  the  resulting  waveform  is  the  BPSK  signal  that  is  easily  demodulated  by  a  conven¬ 
tional  BPSK  receiver.  [20] 

Spread  spectrum  signals  have  some  very  important  benefits.  They  are  very  effi¬ 
cient  at  suppressing  both  multi-user  interference  and  channel-induced  intersymbol  inter¬ 
ference  (ISI)  due  to  multipath  arrivals.  They  can  also  be  used  for  hiding  a  low-power  sig¬ 
nal  below  the  noise  floor.  This  is  important  because  the  detection  of  a  DSSS  signal  by  an 
unauthorized  listener  is  very  difficult.  That  is  the  reason  DSSS  systems  are  referred  to  as 
Low  Probability  of  Detection  (LPD)  communications  systems.  Even  if  this  signal  is  de¬ 
tected,  the  knowledge  of  the  PN  code  is  necessary  in  order  to  demodulate  it.  As  a  result 
DSSS  systems  are  referred  as  Low  Probability  of  Intercept  (LPI)  communications  sys¬ 
tems.  For  the  above  mentioned  reasons,  the  implementation  of  spread  spectrum  in  Sea- 
web  is  desirable.  [21] 

C.  CHANNEL  ESTIMATION 

Consider  a  signal  x(n)  (in  the  discrete  time  domain)  that  is  sent  through  the  un¬ 
derwater  channel.  For  now,  the  format  of  the  signal  is  not  important.  A  usual  way  to 
model  a  multipath  fading  channel  is  to  represent  it  as  a  tapped  delay  line  with  L  taps.  The 
delay  between  the  taps  is  set  equal  to  the  inverse  of  the  sampling  rate,  as 

d=Y~-  <3'2) 

v sample 

The  total  length  L  of  the  tapped  delay  line  identically  has  to  correspond  to  the 
multipath  spread  of  the  channel,  that  is  Tm  =  d  /  L . 

The  received  signal  y{n)  (at  time  n)  is  the  result  of  the  summation  of  the  L  differ¬ 
ent  contributions  of  delayed  versions  of  the  signal  x(n)  weighted  by  the  appropriate 
channel  coefficients.  Then, 

y(n)  =  Yjx(n-k)hk(n)  (3.3) 

k= 0 

The  last  equation  in  vector  form  can  be  written: 
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y(n)  =  h  (n)x(n) 


(3.4) 


where 

x(«) 

x(n-  l) 

x(n-L- 1) 

and 

h(n)  =  (h0(n)  \(n)  •••  hL_x(n)).  (3.6) 

The  function  hk(n)  is  the  impulse  response  of  the  channel,  which  we  introduced 
in  Chapter  II.  This  function  can  be  written  in  polar  form  as  hk(n)  =  ak  (n)ej0t^  where 
ak  (n)  represents  the  amplitude  and  6 k  (n)  the  phase  shift  of  the  underwater  channel  im¬ 
pulse  response  at  time  n  and  time  delay  kd . 

The  method  of  channel  estimation  we  are  going  to  use  in  this  work  uses  a  DSSS 
wideband  signal  in  the  input  signal  x(n) .  In  the  baseband,  this  signal  has  the  form 

x(n)  =  d(n)[pl(n)  +  jp2{n )),  where  px  and p2  are  two  different  PN  Gold  sequences  of 

the  same  length,  modulated  by  the  same  data  bit  (d  (/?)) .  The  period  of  the  PN  sequence 

has  to  be  longer  than  the  multipath  spread  of  the  channel.  Also  the  chip  duration  has  to  be 
such  that  it  will  give  us  a  signal  Bandwidth  BW  =  2jTdup  ,  much  larger  than  the  coherence 

bandwidth  Bc  [Ref.  22],  The  resulting  minimum  time  resolution  in  determining  the  im¬ 
pulse  response  of  the  channel  will  be  Tchip  seconds.  In  this  case  the  received  signal  will  be 
given  by 

L- 1 

r(n)  =  Yjhk  ( n)d (i n)(p ,  ( n-k)  +  jp2  ( n-k )).  (3.7) 

k=0 
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The  receiver  correlates  the  received  signal  with  the  known  spreading  sequence  delayed 
by  kd ,  and  this  process  is  repeated  L  times  for  each  time  delay  from  0  to  (Z-l)  d  sec¬ 
onds.  The  discrete  time  between  transmitted  and  received  sequences  can  be  computed  as: 

n+N-l  1 

c(m,n)=  £  —r(i)(pl(i-m)-jp2(i-m))d(i-m) 

i=n  ^ 

=  ^^-h{m,n){L[p\i-m))  +L(p2{i-m))  }S(i-m )  (3.8) 

1=0  ^ 

N- 1 

=  ^h{rn,n)S(i-m)  =  h(m,n) 

i=0 

where  S(t)  is  the  discrete  time  impulse. 

This  shows  that,  at  least  in  the  ideal  case,  the  impulse  response  of  the  channel  can 
be  computed  by  crosscorrelating  the  transmitted  and  the  received  sequences.  The  assump¬ 
tions  are  that  the  two  sequences,  pl  (n)  and p2  (n) ,  have  ideal  autocorrelations  and  that 

there  is  no  noise  in  the  receiver.  In  practice  the  result  will  be  deteriorated  due  to  the  noise 
and  the  non-ideal  correlation  properties  of  the  Gold  sequence.  [22] 

Using  the  previously  mentioned  method  we  get  an  estimate  of  the  impulse  re¬ 
sponse  of  the  channel  hest  ( m,n ).  It  is  a  function  of  two  variables,  the  time  delay  m  and 

the  absolute  time  n .  By  processing  this  function  we  measure  the  characteristics  of  the 
underwater  channel. 

The  first  function  we  obtain  is  the  multipath  intensity  profile  of  the  channel,  given 
by  the  relation 

(3.9) 

where  Ph  (in)  represents  the  average  power  output  of  the  channel  as  a  function  of  time  de¬ 
lay.  The  width  of  this  function  is  the  multipath  spread  of  the  channel.  As  discussed  in 
Chapter  II  the  coherence  bandwidth  of  the  channel  is  the  inverse  of  the  multipath  spread 
B  =  \/T  . 

c  /  m 
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The  second  function  is  the  Doppler  power  spectrum  of  the  channel,  which  is  the 
Fourier  transform  of  the  spaced- time  correlation  function  of  the  channel 

=  O-IO) 

where  SH  (An)  is  defined  in  Chapter  II.  It  is  used  to  examine  the  Doppler  effects  of  the 

channel,  the  Doppler  spread  and  the  Doppler  shift.  The  Doppler  power  spectrum  Ph  (v) 

is  centered  in  the  frequency  spectrum  in  the  frequency  that  corresponds  to  Doppler  shift 
and  its  bandwidth  corresponds  to  the  Doppler  spread  of  the  underwater  channel.  As  we 
examined  before,  the  coherence  time  of  the  channel  is  the  inverse  of  the  Doppler  spread, 

(H=W 

The  most  useful  function  is  the  third  one,  which  combines  information  about  the 
frequency  and  time  spread  of  the  channel.  It  is  called  the  scattering  function  of  the  chan¬ 
nel.  We  can  get  the  scattering  function  by  taking  the  Fourier  transform  with  respect 
to  An ,  of  the  autocorrelation  function  of  the  estimated  impulse  response 

SH  =  Fu|*(»,i«)),  (3.11) 

where 

1  N 

A  (m’ An)  =  T7 Z  'best  ( m » n)Kst  {rn,n  +  An)}.  (3.12) 

^  n= o 

This  is  a  very  important  function  since  the  knowledge  of  S(m,v )  by  itself  is 

enough  to  give  us  all  the  characteristics  of  the  channel  (Doppler  shift,  Doppler  spread, 
and  multipath  spread).  [23] 

D.  CHAPTER  SUMMARY 

In  this  chapter,  we  developed  our  method  of  channel  estimation.  In  the  next  two 
chapters  we  implement  this  method,  and  apply  it  first  to  an  artificial  channel  (Chapter  IV) 
and  then  to  the  data  of  a  real  ocean  experiment  (Chapter  V). 
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IV.  DEVELOPMENT  OF  METHOD  ON  ARTIFICIAL  CHANNEL 


The  first  objective  in  the  development  and  testing  of  channel  estimation  is  the  ap¬ 
plication  of  the  method  on  an  artificial  time-varying  channel  with  known  characteristics. 
The  scope  of  this  simulation  is  to  confirm  that  the  channel  estimation  algorithm  works 
properly. 

A.  ARTIFICIAL  CHANNEL 

The  received  signal  can  be  represented  as  we  developed  in  the  Chapter  III  by  the 
superposition  of  the  L  delayed  replicas  of  the  transmitted  signal, 

y{t)  =  Yjx{t~k^T)hk{t)  (4-1) 

k= 0 

where  x(t)  is  the  wideband  transmitted  signal  and  hk(t)  is  the  k- th  impulse  response 

component  at  time  t.  The  artificial  channel  is  constructed  according  to  the  structure  of 
Figure  3,  in  which  we  have  three  discrete  paths  for  the  eigenrays.  In  this  representation 
the  transmitter  and  receiver  are  stationary,  which  implies  that  there  is  no  Doppler  shift 
due  to  relative  motion.  Environmental  parameters  such  as  the  roughness  of  the  sea  or  the 
currents  in  the  water  column  produce  a  time- varying  underwater  channel.  Each  path  has  a 
different  time-varying  nature  and  hence  a  different  (nonzero)  Doppler  spread.  The  first 
path  is  the  direct  path,  which  we  assume  to  have  zero  time  delay  and  the  smallest  Doppler 
spread.  The  second  and  third  dominant  paths  correspond  to  the  surface-reflected  eigen¬ 
rays,  having  constant  time  delay  10  and  50  milliseconds,  respectively,  with  different 
time-varying  weights.  The  amplitudes  of  the  three  components  of  the  impulse  response 
are  illustrated  in  Figure  12  as  a  function  of  absolute  time.  As  indicated  in  this  figure, 
there  is  a  substantial  time  variation  in  their  amplitudes  in  a  relatively  small  period  of  time 
(about  4  seconds).  The  average  multipath  intensity  profile  of  the  artificial  channel  as  a 
function  of  the  time  delay  is  illustrated  in  Figure  13;  clearly,  the  multipath  spread  Tm  of 
the  artificial  channel  is  50  milliseconds. 
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E  RESPONSE 


IMPULSE  RESPONSE  AT  TIME  DELAY  ZERO 


Figure  12.  Amplitude  of  the  three  impulse  responses  components  in  absolute  time 


MULTIPATH  INTENSITY  PROFILE  (MIP) 


Figure  13.  Multipath  intensity  profile  of  artificial  signal  (in  time  delay) 
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In  order  to  assess  the  resulting  Doppler  spread  of  the  estimated  channel,  we  de¬ 
velop  the  space-time  correlation  function  of  the  known  artificial  channel,  which  we  intro¬ 
duced  in  Chapter  II.  We  express  this  as  the  autocorrelation  of  each  component  of  the  im¬ 
pulse  response 


Sa(At) 


^{\h*k(t)hk{t  +  M)\2} 
H\K{t)k(t)\2}  '  ’  ’ 


(4.2) 


In  our  case,  the  impulse  response  has  only  three  nonzero  components,  so  we  get  three 
autocorelation  functions,  which  are  illustrated  in  Figure  14.  As  we  discussed  earlier,  the 
period  of  time  over  which  this  function  is  approximately  constant  is  the  coherence  time  of 
the  specific  path.  Let  us  consider  that  the  coherence  time  corresponds  to  the  period  At 
over  which  S*  (At)  is  greater  than  0.95.  The  resulting  coherence  time  for  the  direct  path 
case  is  0.5  seconds,  for  the  second  path  is  0.28  seconds  and  for  the  third  path  is  0.41  sec¬ 
onds.  Since  the  Doppler  spread  Bd  is  equal  to  the  inverse  of  the  coherence  time,  the  cor¬ 
responding  spreads  for  the  three  paths  are  2  Hz,  3.57  Hz  and  2.44  Hz,  respectively. 


COHERENCE  FUNCTION  OF  THE  DIRECT  PATH 


Figure  14.  Coherence  function  of  the  three  components  of  impulse  response 


31 


B.  TRANSMITTED  PROBING  SIGNAL 

The  transmitted  signal  is  a  Direct  Sequence  Spread  Spectrum  (DSSS)  baseband 
signal  like  the  one  studied  in  Chapter  III.  It  uses  PN  sequences  defined  by  two  different 
Gold  codes  of  length  2047,  one  on  the  in-phase  and  the  other  on  the  quadrature  compo¬ 
nent,  modulating  the  same  data  bit.  These  codes  are  taken  from  work  done  by  [24]  and 
the  actual  sequences  were  downloaded  directly  from  the  web  site.  The  chipping  rate  used 
in  the  simulation  is  4000  chips/second.  Each  set  of  2047  chips  represents  one  bit,  so  the 
data  rate  of  the  simulation  will  be  about  2  bps. 

The  length  of  the  Gold  sequence  is  a  very  important  parameter  in  our  analysis.  As 
it  gets  smaller  the  estimation  of  the  channel  characteristics  obviously  degrades.  On  the 
other  hand,  as  it  gets  larger,  the  estimation  of  the  channel  is  more  robust  but,  since  the  re¬ 
sult  of  the  method  is  the  average  of  the  impulse  response  over  the  length  of  the  PN  se¬ 
quence,  it  gets  less  accurate,  so  there  is  a  tradeoff  there.  From  our  experiments,  it  seems 
that  the  choice  of  length  2047  is  a  reasonable  compromise.  In  what  follows,  we  send  an 
information  sequence  of  8  bits,  so  the  duration  of  the  entire  information  sequence  is  4.094 
seconds.  In  the  first  case  we  study  an  ideal  noise-free  case.  In  the  second  case,  the  signal 
is  corrupted  by  the  presence  of  additive  white  Gaussian  noise  (AWGN). 

C.  RESULTS  OF  THE  METHOD 

In  the  following  simulation  the  channel  parameters  are  estimated  without  noise, 
and  later  with  additive  white  Gaussian  noise  in  the  channel. 

1.  Results  for  the  Case  Without  Noise 

In  the  upper  half  of  Figure  15  the  exact  form  of  the  impulse  response  of  the  direct 
path  is  shown,  whereas  in  the  lower  half  is  the  algorithm’s  estimate  for  the  same  function. 
Clearly,  we  do  not  get  the  exact  value  of  the  impulse  response,  but  an  average  estimate  of 
the  next  2047  values  of  the  impulse  response,  which  cover  a  period  of  about  half  a  sec¬ 
ond.  As  a  result  of  that,  it  seems  that  the  estimate  has  a  time  offset  of  250  milliseconds 
from  the  actual  one.  The  same  happens  with  the  other  two  components  of  h(t),  corre¬ 
sponding  to  time  delays  10  and  50  milliseconds,  that  are  illustrated  in  Figures  16  and  17, 
respectively.  Except  for  the  fact  that  they  are  averages  and  not  exact  values,  the  estimates 
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in  all  three  cases  are  quite  accurate,  which  means  that  until  this  stage  the  method  seems  to 
be  accurate. 


ACTUAL  IMPULSE  RESPONSE  OF  THE  DIRECT  PATH 


Figure  15.  Exact  and  estimated  impulse  response  for  the  direct  path 


ACTUAL  IMPULSE  RESPONSE  OF  THE  SECOND  PATH  (TIME  DELAY  10  mecs) 


Figure  16.  Exact  and  estimated  impulse  response  for  the  2nd  path 
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ACTUAL  IMPULSE  RESPONSE  OF  THE  THIRD  PATH  (TIME  DELAY  50  mecs) 


Figure  17.  Exact  and  estimated  impulse  response  for  the  3rd  path 


The  resulting  multipath  intensity  profile  of  the  underwater  channel  is  illustrated  in 
Figure  18.  The  result  is  impressive  since  the  estimate  multipath  intensity  profile  (MIP)  is 
almost  the  same  as  the  original  MIP  illustrated  in  Figure  13. 


ESTIMATED  MULTIPATH  INTENSITY  PROFILE  OF  ARTIFICIAL  UNDERWATER  CHANNEL 


Figure  18.  Estimated  Multipath  intensity  profile  of  artificial  signal  (in  time  delay) 
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The  most  significant  function  is  the  scattering  function  S(t,  v)  of  the  channel 
which  we  analyzed  in  Chapter  III.  We  recall  that  it  combines  information  about  both  the 
frequency  and  time  spread  of  the  channel.  In  Figure  19  the  normalized  scattering  function 
of  the  channel  is  illustrated.  Just  by  inspection  of  the  plot,  we  can  conclude  the  following: 

•  The  channel  has  three  dominant  paths  with  time  delays  0,  10  and  50  milli¬ 
seconds  respectively. 

•  It  seems  that  the  second  path  is  stronger  than  the  other  two,  but  this  is  not 
true.  The  apparent  discrepancy  is  because  we  took  out  the  DC  (zero  Dop¬ 
pler)  component  of  the  impulse  response  before  processing.  We  did  that  for 
presentation  reasons,  so  that  the  Doppler  spread  would  be  more  obvious. 

•  The  Doppler  spread  of  the  second  path  is  by  far  the  largest,  the  next  larger  is 
the  third  path  and  the  smaller  Doppler  spread  corresponds  to  the  direct  path. 

•  The  Doppler  shift  of  the  underwater  channel  is  zero,  since  the  functions  are 
centered  around  the  zero  frequency,  consistent  with  a  fixed  transmitter- 
receiver  geometry. 

ESTIMATED  (NORMALIZED)  SCATTERING  FUNCTION  OF  THE  CHANNEL 


FREQUENCY  IN  Hz 


TIME  DELAY  IN  mseconds 


Figure  19.  Estimated  normalized  scattering  function  of  the  artificial  channel 
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The  Doppler  spread  of  the  dominant  paths  is  estimated  by  finding  the  bandwidths 
of  the  scattering  functions  S(t,  v)  corresponding  to  those  paths.  In  order  to  determine  the 
bandwidth  we  follow  the  definition  [25]: 


rr-\sM)w 

J-  00 _ 

r\s,v)\df 

J-  00 


(4.3) 


The  resulting  estimated  values,  as  well  as  the  values  calculated  earlier  from  the  coherence 
times  (A t)c,  of  Doppler  spreads,  are  summarized  in  Table  1.  We  notice  some  differences 

in  the  values,  which  are  due  to  the  arbitrary  threshold  of  0.95  for  the  coherence  time,  and 
to  the  averaging  which  results  in  a  small  distortion  in  the  impulse  response  and  small  er¬ 
rors  in  the  Doppler  spread. 


Path 

Calculated  Doppler 
spread  from  coherence 
time  of  the  channel 

Estimated  Doppler 
spread  using  the 
algorithm 

Direct  path 

2  Hz 

1.3  Hz 

2nd  path  (10  ms  delay) 

3.57  Hz 

4.3  Hz 

3rd  path  (50  ms  delay) 

2.44  Hz 

2.8  Hz 

Table  1.  Doppler  spreac 

Is  of  the  dominant  paths  of  the  underwater  channel 

2.  Results  for  the  AWGN  Case 

Consider  a  channel  in  which  the  received  signal  is  corrupted  by  additive  white 
Gaussian  noise  as  well.  Using  our  previous  model  and  starting  with  very  low  noise  in  the 
channel  (SNR  of  25  dB),  we  notice  that  there  is  no  any  difference  in  the  estimated  im¬ 
pulse  response.  The  method  seems  also  to  be  robust  at  higher  noise  levels.  So  even  in  15- 
dB  SNR,  the  effect  of  the  noise  is  negligible.  This  effect  becomes  more  observable  (but 
not  disturbing)  at  10-dB  SNR  and  below.  The  robustness  of  the  algorithm  results  from  the 
length  of  the  PN  sequence,  which  in  a  sense,  averages  out  the  noise.  The  plots  in  Figures 
21,  22,  23  refer  to  a  signal-to-noise  ratio  of  7  dB. 
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ACTUAL  IMPULSE  RESPONSE  OF  THE  DIRECT  PATH 


ESTIMATED  IMPULSE  RESPONSE  OF  THE  DIRECT  PATH  (AWGN  CASE  /  SNR=7dB) 


Figure  20.  Exact  and  estimated  impulse  response  for  direct  path  (SNR  =  7  dB) 


ACTUAL  IMPULSE  RESPONSE  OF  THE  SECOND  PATH  (TIME  DELAY  10  mecs) 


ACTUAL  IMPULSE  RESPONSE  OF  THE  SECOND  PATH  (TIME  DELAY  10  mecs)  (AWGN  CASE  /  SNR=7dB) 


Figure  21 .  Exact  and  estimated  impulse  response  for  the  2nd  path  (SNR  =  7  dB) 
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Figure  22.  Exact  and  estimated  impulse  response  for  the  3rd  path  (SNR  =  7  dB) 


Figure  23  shows  the  multipath  intensity  profile  of  the  underwater  noisy  channel. 
The  result  is  impressive,  since  the  plot  is  identical  to  the  case  when  noise  is  not  present. 
Figure  24  illustrates  the  scattering  function  of  the  channel.  The  resulting  Doppler  spread 
values  are  almost  the  same  with  the  noise  free  case  and  they  are  summarized  in  Table  2. 


Path 

Estimated  Doppler 
spread 

Noise  free  case 

Estimated  Doppler  spread 
AWGN  case  (SNR  =  7  dB) 

Direct  path 

1.3  Hz 

1.35  Hz 

2nd  path  (10  ms  delay) 

4.3  Hz 

4.3  Hz 

3rd  path  (50  ms  delay) 

2.8  Hz 

2.9  Hz 

Table  2.  Doppler  spreads  of  the  dominant  paths  of  the  noisy  underwater  channel 
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ESTIMATED  MULTIPATH  INTENSITY  PROFILE  OF  ARTIFICIAL  UNDERWATER  CHANNEL  (AWGN  CASE  /  SNR=7dB) 
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Figure  23.  Estimated  Multipath  intensity  profile  of  artificial  signal 

(AWGN  case/ SNR  =  7  dB) 


ESTIMATED  (NORMALIZED)  SCATTERING  FUNCTION  OF  THE  CHANNEL  (AWGN  CASE  /  SNR=7  dB) 


FREQUENCY  IN  Hz 


TIME  DELAY  IN  mseconds 


Figure  24.  Estimated  normalized  scattering  function  of  the  artificial  noisy  channel 


D.  CHAPTER  SUMMARY 

In  this  chapter  we  created  an  artificial  time-varying  channel  with  known  charac¬ 
teristics.  Then,  by  processing  the  received  signal,  which  had  passed  through  the  channel, 

we  obtained  back  the  characteristics  of  the  artificial  channel.  This  artificial  channel 
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worked  as  a  test  channel  to  verify  that  the  method  works  efficiently.  In  the  next  chapter, 
we  process  data  from  an  actual  experiment  (New  England  Shelf  experiment  -  17-20  April 
2000).  Using  the  same  code,  we  determine  the  characteristics  of  the  channel.  We  describe 
the  conditions  of  this  experiment,  the  format  of  the  data  used  and  also  refer  in  the  prob¬ 
lems  encountered  in  extracting  the  scattering  function  from  the  set  of  real  data. 
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V.  NEW  ENGLAND  SHELF  CHANNEL  ESTIMATION 


In  this  chapter  we  validate  the  results  of  this  research  by  using  experimental  data. 
We  measure  the  underwater  channel  parameters  described  theoretically  in  the  Chapter  II. 
Our  measurements  follow  the  method  developed  in  the  Chapter  III  and  tested  for  accu¬ 
racy  in  Chapter  IV. 

A.  DESCRIPTION  OF  THE  NEW  ENGLAND  SHELF  EXPERIMENT 

During  April  2000,  the  experiment  ForeFRONT-2  was  conducted  over  the  New 
England  Shelf  area  charted  in  Figure  25  [26],  As  part  of  this  experiment  the  Signalex-B 
event  obtained  channel  data  [27,  28]. 


Figure  25.  Overview  of  the  Forefront-2  experiment  site. 
The  50-meter  isobath  is  plotted 
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As  illustrated  in  Figure  26,  the  receiver  (telesonar  testbed  on  the  right  side  of  the 
figure  [29])  was  deployed  on  the  ocean  bottom  at  a  depth  of  about  30  meters  in  a  station¬ 
ary  position,  and  recorded  waveforms  (probe  signals)  were  transmitted  from  the  R/V 
Connecticut  as  it  drifted  west  southwest  from  the  receiver.  The  transmitter  was  an  over- 
the-side  projector  deployed  at  a  depth  of  about  20  meters.  The  received  waveforms  were 
taken  at  ten  transmitter-to-receiver  ranges  starting  from  700  meters  and  increasing  to 
6550  meters.  The  purpose  of  these  multiple  ranges  was  to  examine  in  more  detail  the  be¬ 
havior  of  the  underwater  channel  as  the  distance  grows. 


Figure  26.  Rough  illustration  of  the  experiment 


A  plot  of  the  bottom  topography  and  source  track  is  provided  in  Figure  27.  The 
position  of  the  stationary  receiver  is  indicated  with  the  letter  R.  The  transmitter  sent  the 
waveforms  at  10  different  distances  from  the  test  bed.  Those  positions  are  indicated  in 
Figure  27  with  the  letters  T1  to  T10. 
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Figure  27.  Transmitter-  Receiver  positions  during  the  experiment 


In  order  to  determine  the  sound  speed  profile  in  the  underwater  channel,  we  used 
measurements  from  a  CTD  instrument,  where  CTD  stands  for  Conductivity,  Temperature 
and  Depth.  The  measurement  took  place  at  time  20:00,  whereas  the  waveforms  were  re¬ 
ceived  from  21:40  until  02:10,  so  there  is  a  great  possibility  that  the  sound  speed  profile 
changed  slightly  during  this  period  of  time.  All  the  eigenrays  estimation  plots  and  the 
theoretical  impulse  response  results  are  based  on  this  measured  sound  speed  profile, 
which  is  illustrated  in  Figure  28. 
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SOUND  SPEED  PROFILE 


Figure  28.  Sound  speed  profile 


B.  DESCRIPTION  OF  THE  PROBE  SINGAL 

During  the  experiment,  the  transmitter  generated  a  special  type  of  waveform.  The 
waveform  sent  is  shown  in  Figure  29  and  consists  of  different  types  of  probes  [30,  31]. 
The  total  duration  of  the  test  probes  was  about  2  minutes  and  5  seconds.  The  analytic  de¬ 
scription  of  the  various  probes  follows. 


Figure  29.  Transmitted  probe  waveform 
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The  first  probe  consists  of  40  LFM  chirps  with  duration  of  50  milliseconds  each, 
with  200  milliseconds  silence  separating  adjacent  chips  for  a  total  duration  of  10  seconds. 
Each  chirp  sweeps  through  frequencies  from  8  to  16  kHz.  As  we  will  see  later,  this  type 
of  probe  is  useful  in  determining  the  impulse  response  accurately,  but  it  does  not  yield  a 
good  measurement  of  the  scattering  function  (specifically  Doppler  spread)  of  the  channel. 
After  one  second  of  silence,  the  second  probe  is  sent.  It  is  a  comb  of  17  tones  of  20  sec¬ 
onds  duration.  The  17  tones  have  500  Hz  of  separation;  the  lowest  is  positioned  at  8  kHz 
and  the  highest  at  16  kHz.  They  are  useful  in  measuring  the  Doppler  shift  of  the  underwa¬ 
ter  channel.  However,  the  scattering  function  cannot  be  derived  because  as  we  discussed 
earlier,  each  path  has  its  own  Doppler  shift,  so  the  result  we  get  from  the  tone  combina¬ 
tion  is  just  a  rough  estimate  of  the  average  Doppler  shift. 

The  next  two  probe  waveforms  were  not  used  in  the  channel  estimation  in  this 
work.  In  both  cases,  m-sequences  are  transmitted.  In  the  first  case,  ten  m-sequences  are 
sent,  with  m  =  6  and  silent  gaps  in  between  them  are  sent.  In  the  second  case,  ten  m- 
sequences  with  m=  10  are  sent  but  with  no  silent  gaps  in  between  them.  The  chip  rate  is 
4000  chips/second  and  the  center  frequency  is  12  kHz.  The  next  probe,  which  also  was 
not  studied  in  this  work,  follows  after  a  silence  gap  of  1.5  second.  It  is  a  music  clip  band- 
shifted  to  a  center  frequency  of  12  kHz.  The  next  probe  is  the  one  we  are  most  interested 
in.  It  has  almost  the  same  format  with  the  DSSS  signal  we  discussed  in  Chapter  IV.  This 
signal  is  a  bandpass  direct-sequence  spread  spectrum  signal  with  center  frequency  of  12 
kHz,  a  chip  rate  of  2000  chips/second  and  a  sampling  rate  of  48000  samples/second.  Four 
hundred  information  bits  were  transmitted.  We  used  two  different  Gold  sequences  for 
spreading  sequences;  a  different  PN  sequence  is  used  for  the  in-phase  and  quadrature 
data,  although  the  same  information  bit  is  modulated  by  each  component.  Following  Fig¬ 
ure  29,  we  notice  that  we  send  six  different  DSSS  probes.  The  six  different  probes  corre¬ 
spond  to  six  different  bit  rates  of  10,  20,  50,  100,  200  and  400  bps.  In  all  cases,  we  keep 
the  chip  rate  and  length  of  Gold  sequence  constant.  This  means  that,  for  example,  in  the 
first  case  about  400  chips  are  modulated  by  one  data  bit,  in  the  second  case  200  chips  are 
modulated  by  one  data  bit,  etc.  Intuitively,  this  implementation  would  result  in  greater 
correlation  noise  than  the  ideal  case  we  considered  in  Chapter  IV  in  which  the  entire 
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length  of  the  Gold  sequence  was  modulated  by  one  data  bit.  In  order  to  minimize  this 
problem,  we  use  the  lower  bit  rate  probe  (first  probe  with  10  bps)  and  for  the  channel  es¬ 
timation,  we  will  use  a  portion  of  the  Gold  code  and  not  the  entire  sequence.  In  Figure  30, 
the  most  significant  probes  are  summarized  by  illustrating  their  frequency  spectra. 


Figure  30.  Spectrograms  of  the  probes 


C.  CHANNEL  ESTIMATION  RESULTS 

We  now  present  the  results  of  the  channel  estimation  method  for  ten  different 
ranges  and  geometries.  The  first  reception  was  at  time  2143  and  distance  700  meters. 
Subsequent  receptions  occur  every  30  minutes  at  increasing  ranges  until  the  last  one  at 
time  0213  and  at  a  distance  of  6550  meters.  The  ten  cases  are  summarized  in  Table  3. 
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CASE 

TIME 

RANGE  (meters) 

SNR  (dB) 

1 

21:43 

700 

37 

2 

22:13 

1100 

35 

3 

22:43 

1650 

31 

4 

23:13 

2300 

30 

5 

23:43 

3050 

29.5 

6 

00:13 

3700 

25 

7 

00:43 

4350 

25 

8 

01:13 

5000 

16 

9 

01:43 

5700 

13 

10 

02:13 

6550 

7 

Table  3.  Summary  of  the  10  cases 


1.  Received  Signal  at  a  Distance  of  700  Meters 

The  signal  received  at  time  2143  and  at  a  distance  of  700  meters  is  plotted  in  Fig¬ 
ure  3 1 .  At  this  distance,  the  signal  is  very  strong  and  the  resulting  approximate  signal-to- 
noise  ratio  is  37  dB.  The  estimation  of  the  signal-to-noise  ratio  is  done  by  determining 
first  the  noise  power  level  (during  a  time  period  when  the  signal  is  not  present),  and  then 
the  signal  plus  noise  power  level,  during  the  time  period  when  the  DSSS  signal  is  on.  By 
processing  those  two  values,  we  determine  the  SNR  using  the  relation: 


SNRJll^zl 

N 


(5.1) 


Figure  3 1 .  Received  signal  at  distance  of  700  meters 
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In  Figure  32  the  eigenrays  characterizing  the  propagation  channel  are  traced  using 
a  numerical  model  called  Bellhop  [32],  It  is  interesting  to  notice  that  the  first  two  arrivals 
are  almost  simultaneous.  The  other  6  eigenrays  come  into  the  receiver  in  pairs,  and  this  is 
due  to  the  proximity  of  the  receiver  to  the  seabed. 


BELLHOP-  ENV  Ol 
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Range  (m) 

Figure  32.  Eigenrays  plot  for  distance  of  700  meters 

We  first  estimate  the  impulse  response  of  the  channel  by  applying  a  matched  filter 
to  the  probe  pulses  (using  the  known  50  milliseconds  LFM  chirp  as  the  correlation  ker¬ 
nel).  This  method  produces  accurate  measurements  of  the  channel  impulse  response  at 
the  pulse  repetition  rate  (i.e,  4  measurements/second).  The  results  for  the  distance  of  700 
meters  are  shown  in  Figure  33.  We  next  use  the  DSSS  signal  in  the  bandpass,  following 
the  method  we  discussed  in  Chapters  III  and  IV.  The  results  are  shown  in  Figure  34.  Fi¬ 
nally,  we  determine  the  impulse  response  of  the  channel  using  the  DSSS  in  the  baseband. 
This  method  involves  the  complicated  process  of  shifting  the  waveform  from  bandpass  to 
baseband.  The  result  is  shown  in  Figure  35.  The  last  two  methods  give  us  a  measurement 
rate  of  about  24  Hz. 
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Comparing  the  three  plots  of  impulse  response  estimates  we  can  see  that  all  meth¬ 
ods  give  the  same  estimate,  but  for  DSSS,  the  noise  level  is  much  higher  than  when  using 
LFM.  The  explanation  for  this  is  the  high  correlation  noise  we  discussed  previously.  On 
the  other  hand,  the  low  measurement  rate  in  the  LFM  case  will  not  allow  us  to  derive  the 
scattering  function  of  the  channel,  since  aliasing  due  to  low  measurement  rate  does  not 
capture  the  channel  variability.  In  the  other  two  cases,  the  sampling  rate  of  24  Hz  is  fast 
enough  for  this  channel. 

The  last  issue  we  need  to  comment  on  is  the  slope  of  the  estimated  impulse  re¬ 
sponse  shown  in  Figures  33,  34  and  35.  This  is  the  result  of  the  drifting  of  the 
boat/transmitter  and  it  is  an  indication  of  the  Doppler  shift  in  the  signal.  The  greater  the 
slope,  the  larger  the  Doppler  shift. 


|h|2  USE  OF  LFM  chirps 


6^ 


Figure  33.  Estimated  impulse  response  at  700  meters  using  LFM  chirp 
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|h|2  USE  OF  DSSS  SIGNAL  IN  BANDPASS 


x  10-3 


Figure  34.  Estimated  impulse  response  at  700  meters  using  DSSS  signal  in  bandpass 


Figure  35.  Estimated  impulse  response  at  700  meters  using  DSSS  signal  in  baseband 

Using  the  DSSS  signal  at  baseband,  we  now  derive  the  average  multipath  inten¬ 
sity  profile  of  the  channel.  A  preliminary  step  is  to  align  the  impulses  to  compensate  for 
motion.  There  are  many  techniques  to  do  this,  but  the  one  we  used  is  probably  the  most 
accurate  and  fastest.  Starting  from  the  first  estimate,  we  find  the  first  highest  peak;  then 
we  look  in  the  next  sample  to  find  the  highest  peak  searching  only  in  the  neighborhood  of 
the  peak  of  the  previous  sample.  The  method  worked  precisely  for  all  cases.  The  multi- 
path  intensity  profile  Ph  ( r)  was  analyzed  in  Chapter  II.  The  resulting  multipath  intesity 

profile  for  a  distance  of  700  meters  is  illustrated  in  Figure  36.  From  this  figure,  we  can 

determine  the  multipath  spread  of  the  channel  which  is  about  18  milliseconds.  The  Bell- 
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hop  numerical  method  estimate  of  the  MIP  is  presented  in  Figure  37.  Note  that  our  Bell¬ 
hop  channel  modeling  neglects  boundary  losses  and  shows  arrivals  with  artificially  high 
intensity.  The  matching  of  the  measurement  with  the  theoretical  results  is  satisfactory  in 
terms  of  multipath  arrival  structure. 


Figure  36.  Multipath  Intensity  Profile  at  700  meters  using  DSSS  signal  in  baseband 


THEORE^If|AL  MULTIPATH  hNTENSITY'  PROFIT  USING  BbLLHUP  PROGRAM  AT  DISTANCE  7QQ  METERS 


Figure  37.  Multipath  Intensity  Profile  at  700  meters  -  Bellhop  theoretical  estimate 

Following  the  method  of  Chapter  III  using  the  DSSS  signal  in  baseband,  we  de¬ 
termine  the  scattering  function  of  the  underwater  channel  (shown  in  Figure  38)  for  a  dis¬ 
tance  of  700  meters.  The  impressive  and  expected  result  is  that  from  a  noisy  impulse  re¬ 
sponse,  we  get  an  unambiguous  scattering  function  plot.  Clearly,  there  is  a  negative  Dop¬ 
pler  shift  consistent  with  the  opening  range  of  the  source-receiver  geometry  and  a  Dop¬ 
pler  spread  different  for  each  path.  Following  Figure  36,  we  can  see  that  there  are  nine 
dominant  paths  (those  with  magnitude  greater  than  5%  of  the  magnitude  of  the  stronger 
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path).  The  estimated  values  of  Doppler  shift  and  Doppler  spread  for  those  dominant  paths 
are  shown  in  Table  4.  For  comparison  purposes,  we  used  the  tone  probes  to  determine  the 
Doppler  shift,  and  we  got  a  value  of  2.1  Hz.  This  is  quite  close  to  1.9  Hz  which  is  the 
value  of  the  Doppler  shift  corresponding  to  the  first  path.  Another  interesting  result  is  that 
the  stronger  paths  have  relatively  smaller  Doppler  spreads  of  1.3  to  1.6  Hz  when  the 
weaker  paths  have  larger  spreads  on  the  order  of  2.8  Hz.  This  makes  sense  because  the 
weaker  paths  are  experience  more  reflections  at  the  nonstationary  sea  surface  than  the 
stronger  paths,  so  we  should  expect  them  to  have  a  higher  Doppler  spread.  Also,  accord¬ 
ing  to  some  researchers,  this  high  variability  of  the  weaker  paths  is  what  can  cause  an 
equalizer  (for  example  a  Feedback  Decision  Equalizer  using  Recursive  Least  Squares  al¬ 
gorithm)  to  fail  in  a  coherent  reception  communication  scheme.  [33,  34,  35,  36] 


SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  700  METERS 


Figure  38.  Estimated  Scattering  function  of  the  channel  at  distance  of  700  meters 
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Table  4. 


PATH  WITH  TIME 

DOPPLER 

DOPPLER 

DELAY  (in  msec): 

SHIFT  (HZ) 

SPREAD  (HZ) 

0 

-1.9 

1.27 

2.25 

-1.9 

1.50 

2.75 

-1.85 

1.53 

4.375 

-1.75 

1.76 

5.125 

-1.6 

2.30 

11.875 

-1.65 

2.19 

12.875 

-1.65 

2.02 

16.125 

-1.2 

2.84 

17.25 

-1.4 

2.8 

Doppler  spreads  and  shifts  of  the  dominant  paths  at  distance  of  700  meters 


2.  Received  Signal  at  a  Distance  of  1100  Meters 

The  signal  received  at  time  2213  and  at  a  distance  of  1 100  meters  is  illustrated  in 
Figure  39.  As  we  can  see,  the  signal  is  very  strong  at  this  distance  as  well,  and  the  result¬ 
ing  approximate  signal-to-noise  ratio  comes  out  to  be  35  dB.  The  eigenrays’  propagation 
is  traced  in  Figure  40.  We  see  again  that  the  eigenrays  come  into  the  receiver  in  pairs. 
The  estimated  impulse  responses  using  the  three  different  methods  are  illustrated  in  Fig¬ 
ures  41,  42  and  43.  The  results  in  all  three  cases  are  comparable.  The  slope  in  this  case  is 
greater  than  before,  which  means  that  the  boat  is  drifting  faster  and  the  resulting  Doppler 
shift  is  larger. 


Figure  39.  Received  signal  at  distance  of  1 100  meters 
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Figure  40.  Eigenrays  plot  for  distance  of  1 100  meters 


|h|2  USE  OF  LFM  chirps  at  2213 


Figure  41 .  Estimated  impulse  response  at  1100  meters 

using  LFM  chirp 


|h|2  USE  OF  DSSS  SIGNAL  IN  BANDPASS  x104 


Figure  42.  Estimated  impulse  response  at  1100  meters 
using  DSSS  signal  in  bandpass 
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|h|2  USE  OF  DSSS  SIGNAL  IN  BASEBAND 


Figure  43.  Estimated  impulse  response  at  1100  meters 
using  DSSS  signal  in  baseband 


Figure  44  shows  the  multipath  intensity  profile  at  this  range.  We  can  distinguish 
eight  dominant  paths  in  the  underwater  channel.  Also,  the  multipath  spread  is  smaller 
than  before,  about  12  milliseconds.  The  multipath  intensity  profile  as  a  result  of  the  Bell¬ 
hop  model  is  illustrated  in  Figure  45.  The  matching  between  theory  and  observations  is 
still  satisfactory.  In  Figure  46,  the  clear  plot  of  the  scattering  function  is  shown.  Just  by 
inspection,  we  observe  that  in  this  case  the  Doppler  shift  is  larger  and  the  Doppler  spread 
is  smaller  than  before.  Indeed  in  Table  5,  where  the  resulting  Doppler  spreads  and  shifts 
are  presented  for  the  dominant  paths,  we  can  see  that,  for  the  first  path,  the  Doppler  shift 
is  negative  2.25  Hz  (larger  than  before)  and  the  Doppler  spread  is  0.83  Hz  (smaller  than 
in  the  case  at  700  meters).  The  tone  probe  result  for  the  Doppler  shift  is  2.3  Hz  as  well. 
The  maximum  Doppler  spread  has  a  value  of  2.62  Hz. 
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Figure  44.  Multipath  Intensity  Profile  at  1100  meters  using  DSSS  signal  in  baseband 


Figure  45.  Multipath  Intensity  Profile  at  1100  meters  -  Bellhop  theoretical  estimate 


SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  1100  METERS 
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Figure  46.  Estimated  Scattering  function  of  the  channel  at  distance  of  1 100  meters 


56 


PATH  WITH  TIME 

DOPPLER 

DOPPLER 

DELAY  (in  msec): 

SHIFT  (HZ) 

SPREAD  (HZ) 

0 

-2.25 

0.83 

1.5 

-2.15 

1.50 

2.875 

-2.05 

1.70 

3.25 

-2.1 

1.84 

7.625 

-1.8 

2.25 

8.125 

-1.8 

2.13 

10.75 

-1.65 

2.62 

11.375 

-1.6 

2.40 

Table  5.  Doppler  spreads  and  shifts  of  the  dominant  paths  for  distance  of  1 100  meters 

3.  Received  Signal  at  a  Distance  of  1650  Meters 

The  signal  received  at  time  2243  and  at  a  distance  of  1650  meters  is  illustrated  in 
Figure  47.  As  we  can  see,  the  signal  is  strong  at  this  distance  as  well,  and  the  resulting 
signal-to-noise  ratio  is  approximately  31  dB.  In  Figure  48  the  eigenrays’  propagation, 
modeled  by  the  Bellhop  numerical  analysis,  is  traced.  In  this  case  too,  the  rays  arrive  in 
pairs  at  the  receiver,  but  the  effect  of  attenuation  is  much  stronger,  and  in  the  actual  im¬ 
pulse  response  plot  those  pairs  are  less  distinguishable.  The  estimated  impulse  responses 
using  the  three  different  methods  are  presented  in  Figures  49,  50  and  51.  The  results  in  all 
three  cases  are  consistent.  The  slope  in  this  case  is  even  greater  than  in  two  previous 
cases,  which  means  that  the  boat  is  drifting  even  faster  and  the  resulting  Doppler  shift  is 
larger. 
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Figure  47.  Received  signal  at  distance  of  1650  meters 


Figure  48.  Eigenrays  plot  for  distance  of  1650  meters 


Figure  49.  Estimated  impulse  response  at  1650  meters 

using  LFM  chirp 


Figure  50.  Estimated  impulse  response  at  1650  meters 
using  DSSS  signal  in  bandpass 


Figure  5 1 .  Estimated  impulse  response  at  1650  meters 


using  DSSS  signal  in  baseband 
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Figure  52  shows  the  multipath  intensity  profile  at  1650  meters.  At  this  distance, 
we  can  distinguish  six  dominant  paths  in  the  underwater  channel  and  the  multipath  spread 
is  approximately  13  milliseconds.  The  multipath  intensity  profile  as  a  result  of  the  Bell¬ 
hop  program  is  given  in  Figure  53.  As  we  can  see,  the  attenuation  effects  are  not  captured 
by  the  numerical  simulation.  It  just  gives  us  a  feeling  about  the  multipath  delays.  Figure 
54  shows  the  plot  of  the  scattering  function  of  the  channel.  Just  by  inspection,  we  can  ob¬ 
serve  that  in  this  case  the  Doppler  shift  is  larger  than  in  the  two  previous  cases;  also  the 
Doppler  spread  is  quite  large  as  well.  In  Table  6,  the  estimated  Doppler  spreads  and  shifts 
are  shown  for  the  dominant  paths.  The  estimated  Doppler  spread  varies  from  1.38  Hz 
(stronger  path)  to  3.35  Hz  (weaker  path).  The  estimated  Doppler  shift  of  the  first  path  is 
indeed  larger  than  before  and  it  has  a  value  of  negative  3.4  Hz.  By  comparison,  the  tone 
combination  method  results  in  a  negative  Doppler  shift  of  3.1  Hz. 


Figure  52.  Multipath  Intensity  Profile  at  1650  meters  using  DSSS  signal  in  baseband 


THEORE*l$0kL  MULTIPATH  INTENSITY  PROFILE  USING  BELLHOP  PROGRAM  AT  DISTANCE  1650  METERS 
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Figure  53.  Multipath  Intensity  Profile  at  1650  meters  -  Bellhop  theoretical  estimate 
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SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  1650  METERS 


Figure  54.  Estimated  Scattering  function  of  the  channel  at  distance  of  1650  meters 


PATH  WITH  TIME 

DOPPLER 

DOPPLER 

DELAY  (in  msec): 

SHIFT  (HZ) 

SPREAD  (HZ) 

0 

-3.4 

1.38 

2 

-3.4 

2.24 

4.875 

-2.7 

2.76 

5.375 

-3 

2.2 

8.875 

-2.5 

2.98 

13 

-2 

3.35 

Table  6.  Doppler  spreads  and  shifts  of  the  dominant  paths  for  distance  of  1650  meters 


4.  Received  Signal  at  a  Distance  of  2300  Meters 

The  signal  received  at  time  2313  and  at  a  distance  of  2300  meters  is  plotted  in 
Figure  55.  Looking  at  the  waveform,  we  notice  that  the  signal  is  still  quite  strong  result¬ 
ing  in  a  signal-to-noise  ratio  of  30  dB.  In  Figure  56  the  eigenray  propagation  for  the  dis¬ 
tance  of  2300  meters  is  modeled  and  traced.  In  this  case  the  propagation  becomes  very 
complicated  due  to  the  increased  distance.  The  estimated  impulse  responses  using  the 
three  different  methods  are  presented  in  Figures  57,  58  and  59. 
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Figure  55.  Received  signal  at  distance  of  2300  meters 


BELLHOP-  ENV  01 


Figure  56.  Eigenrays  plot  for  distance  of  2300  meters 
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Figure  57.  Estimated  impulse  response  at  2300  meters 

using  LFM  chirp 


Figure  58.  Estimated  impulse  response  at  2300  meters 
using  DSSS  signal  in  bandpass 


Figure  59.  Estimated  impulse  response  at  2300  meters 
using  DSSS  signal  in  baseband 
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The  multipath  intensity  profile  at  2300  meters  is  shown  in  Figure  60.  At  this  dis¬ 
tance,  we  can  distinguish  eight  dominant  paths  in  the  underwater  channel,  and  the  multi- 
path  spread  is  about  10.5  milliseconds.  The  multipath  intensity  profile  as  a  result  of  the 
Bellhop  model  is  illustrated  in  Figure  61.  The  result  is  quite  similar  except  that  the  paths 
corresponding  to  time  delays  3.25  and  3.625  msecs  are  absent  from  the  model.  Figure  62 
indicates  the  scattering  function  of  the  channel.  Just  by  inspection  we  conclude  that  there 
is  a  strong  negative  Doppler  shift  and  a  few  dominant  paths.  In  Table  7,  the  estimated 
Doppler  spreads  and  shifts  are  presented  for  the  dominant  paths.  The  estimated  Doppler 
spread  varies  from  0.89  Hz  (stronger  path)  to  2.85  Hz  (weaker  path).  We  estimate  that  the 
first  path  has  a  negative  Doppler  shift  with  a  value  of  3.2  Hz.  By  comparison,  the  tone 
combination  method  results  in  a  negative  Doppler  shift  of  3.1  Hz. 


Figure  60.  Multipath  Intensity  Profile  at  2300  meters  using  DSSS  signal  in  baseband 


THEOREtlfl»t  MULTIPATH  INTENSITY  PROFILE  USING  BELLHOP  PROGRAM  AT  DISTANCE  2300  METERS 

4 

3.5 

3 


2.5 

i 

1 


TIME  DELAY  IN  MILLISECONDS 


Figure  61.  Multipath  Intensity  Profile  at  2300  meters  -  Bellhop  theoretical  estimate 
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SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  2300  METERS 
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Figure  62.  Estimated  Scattering  function  of  the  channel  at  distance  of  2300  meters 


Table  7. 


PATH  WITH  TIME 

DOPPLER 

DOPPLER 

DELAY  (in  msec): 

SHIFT  (HZ) 

SPREAD  (HZ) 

0 

-3.2 

0.89 

0.75 

-3.2 

1.10 

1.5 

-2.9 

2.00 

3.25 

-2.5 

2.56 

3.625 

-2.4 

2.85 

5.625 

-2.6 

2.45 

6 

-2.6 

2.40 

10.375 

-2.5 

2.80 

Doppler  spreads  and  shifts  of  the  dominant  paths  for  distance  of  2300  meters 


5.  Received  Signal  at  a  Distance  of  3050  Meters 

The  signal  received  at  time  2343  and  at  a  distance  of  3050  meters  is  plotted  in 
Figure  63.  The  signal  to  noise  ratio  is  not  very  different  than  before  with  a  value  of  29.5 
dB.  In  Figure  64  the  eigenray  propagation  for  the  distance  of  3050  meters,  a  result  of  the 
Bellhop  numerical  analysis,  is  traced.  There  is  still  more  than  one  direct  path  and  many 
reflected  paths  in  the  propagation.  The  estimated  impulse  response  functions  using  the 
three  different  methods  are  presented  in  Figures  65,  66  and  67. 
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Figure  63.  Received  signal  at  distance  of  3050  meters 
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Figure  64.  Eigenrays  plot  for  distance  of  3050  meters 
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Figure  65.  Estimated  impulse  response  at  3050  meters 

using  LFM  chirp 


Figure  66.  Estimated  impulse  response  at  3050  meters 
using  DSSS  signal  in  bandpass 


Figure  67.  Estimated  impulse  response  at  3050  meters 
using  DSSS  signal  in  baseband 
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The  multipath  intensity  profile  at  3050  meters  is  plotted  in  Figure  68.  It  is  inter¬ 
esting  to  notice  that  at  this  distance,  the  multipath  effect  is  almost  negligible  since  there 
are  only  two  dominant  paths  in  the  underwater  channel.  More  specifically,  the  channel 
has  one  strong  path  and  one  very  weak  (its  intensity  is  equal  to  one  tenth  of  that  of  the 
strong),  so  the  multipath  spread  is  just  0.6  milliseconds.  The  Bellhop  model  shown  in 
Figure  69  shows  many  arrivals  after  the  measured  multipath  spread  of  the  channel.  Proba¬ 
bly  those  arrivals  are  very  weak  due  to  the  multiple  reflections.  Figure  70  indicates  the 
scattering  function  of  the  channel  with  the  two  paths.  In  Table  8,  the  estimated  Doppler 
spreads  and  shifts  are  presented  for  the  two  paths.  The  estimated  Doppler  spreads  are 
0.53  and  1.44  corresponding  to  the  strong  and  weak  path  of  the  channel,  respectively.  The 
Doppler  shift  of  the  first  path  is  negative  with  a  value  of  3.5  Hz.  By  comparison,  the  tone 
combination  method  results  in  a  negative  Doppler  shift  of  3.3  Hz. 


1  O  '■  MULTIPATH  INTENSITY  PROFILE  AT  DISTANCE  3050  METERS 


Figure  68.  Multipath  Intensity  Profde  at  3050  meters  using  DSSS  signal  in  baseband 


Figure  69.  Multipath  Intensity  Profile  at  3050  meters  -  Bellhop  theoretical  estimate 
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SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  3050  METERS 
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Figure  70. 


Estimated  Scattering  function  of  the  channel  at  distance  of  3050  meters 


Table  8. 


PATH  WITH  TIME 
DELAY  (in  msec): 

DOPPLER 
SHIFT  (HZ) 

DOPPLER 
SPREAD  (HZ) 

0 

-3.5 

0.53 

0.625 

-3.2 

1.44 

Doppler  spreads  and  shi 


'ts  of  the  dominant  paths  for  distance  of  3050  meters 


6.  Received  Signal  at  a  Distance  of  3700  Meters 

The  signal  received  at  time  0013  and  at  a  distance  of  3700  meters  is  presented  in 
Figure  71.  The  signal-to-noise  ratio  is  25  dB.  In  Figure  72  the  eigenray  propagation  for 
the  distance  of  3700  meters  predicted  by  the  Bellhop  numerical  analysis,  is  traced.  The 
complexity  of  the  propagation  in  the  distance  of  3700  meters  is  very  high.  The  estimated 
impulse  response  functions  using  the  three  different  methods  are  presented  in  Figures  73, 
74  and  75. 
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Figure  7 1 .  Received  signal  at  distance  of  3700  meters 
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Figure  72.  Eigenrays  plot  for  distance  of  3700  meters 
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Figure  73.  Estimated  impulse  response  at  3700  meters 

using  LFM  chirp 


Figure  74.  Estimated  impulse  response  at  3700  meters 
using  DSSS  signal  in  bandpass 


Figure  75.  Estimated  impulse  response  at  3700  meters 
using  DSSS  signal  in  baseband 
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The  multipath  intensity  profile  at  3700  meters  is  illustrated  in  Figure  76.  The  di¬ 
rect  path  is  weaker,  and  the  relative  strengths  of  the  later  paths  are  more  significant.  We 
have  four  dominant  paths  (two  strong  and  two  weak).  The  resulting  multipath  spread  has 
a  value  of  3.5  milliseconds.  The  Bellhop  model  shown  in  Figure  77  has  a  good  match  in 
the  first  4  milliseconds  but,  after  that  value  of  delay,  it  shows  3  arrivals  that  we  cannot 
detect  in  the  measured  result.  It  seems  that  they  attenuated  due  to  the  multiple  reflections. 
In  Figure  78,  the  scattering  function  for  this  case  is  illustrated.  In  Table  9,  the  estimated 
Doppler  spreads  and  shifts  are  presented  for  the  four  paths.  The  estimated  Doppler 
spreads  vary  from  0.93  to  2.71  Hz,  corresponding  to  the  strongest  and  the  weakest  path, 
respectively.  The  Doppler  shift  of  the  first  path  is  negative  with  a  value  of  3.5  Hz.  By 
comparison,  the  tone  probes  show  a  negative  Doppler  shift  of  3.6  Hz. 


Figure  76.  Multipath  Intensity  Profile  at  3700  meters  using  DSSS  signal  in  baseband 


Figure  77.  Multipath  Intensity  Profile  at  3700  meters  -  Bellhop  theoretical  estimate 
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SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  3700  METERS 


Figure  78.  Estimated  Scattering  function  of  the  channel  at  distance  of  3700  meters 


PATH  WITH  TIME 

DOPPLER 

DOPPLER 

DELAY  (in  msec): 

SHIFT  (HZ) 

SPREAD  (HZ) 

0 

-3.5 

0.93 

1.5 

-3.3 

1.43 

3.125 

-2.5 

2.71 

3.5 

-2.8 

2.55 

Table  9.  Doppler  spreads  and  shifts  of  the  dominant  paths  for  distance  of  3700  meters 


7.  Received  Signal  at  a  Distance  of  4350  Meters 

The  signal  received  at  time  0043  and  at  a  distance  of  4350  meters  is  illustrated  in 
Figure  79.  Even  through  the  distance  increased  by  650  meters  from  the  last  measurement, 
the  signal-to-noise  ratio  has  the  same  value  as  previously,  25  dB.  The  eigenray  propaga¬ 
tion  modeled  by  the  Bellhop  program  for  the  distance  of  4350  meters  is  presented  in  Fig¬ 
ure  80.  There  is  still  a  direct  path  in  the  propagation  due  to  the  weak  ducting  of  the  eigen- 
rays  in  the  channel. 
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Figure  79.  Received  signal  at  distance  of  4350  meters 
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Figure  80.  Eigenrays  plot  for  distance  of  4350  meters 
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Figure  8 1 .  Estimated  impulse  response  at  4350  meters 

using  LFM  chirp 
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Figure  82.  Estimated  impulse  response  at  4350  meters  using  DSSS  signal  in  bandpass 


Figure  83.  Estimated  impulse  response  at  4350  meters  using  DSSS  signal  in  baseband 
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The  multipath  intensity  profile  of  the  channel  at  4350  meters  is  illustrated  in  Fig¬ 
ure  84.  The  multipath  effect  at  this  distance  is  getting  very  weak  again  since,  even  if  we 
have  three  dominant  paths  in  this  case,  the  relative  intensities  of  the  second  and  third 
paths,  in  contrast  with  the  first  path,  are  very  small.  The  resulting  multipath  spread  has  a 
value  of  2.4  milliseconds.  The  Bellhop  model  shown  in  Figure  85  does  a  very  good  fit  to 
the  measured  result,  except  the  peak  at  8.7  msecs  delay  which  is  probably  attenuated  in 
the  real  propagation.  In  Figure  86,  the  scattering  function  of  the  channel  is  illustrated.  In 
Table  10,  we  summarize  the  estimated  Doppler  spreads  and  shifts  of  the  three  paths.  The 
estimated  Doppler  spreads  vary  from  0.66  to  2.06  Hz,  corresponding  to  the  strongest  and 
the  weakest  path  respectively.  The  Doppler  shift  of  the  first  path  is  negative  with  a  value 
of  3. 1 5  Hz.  For  comparison,  the  tone  probe  show  a  negative  Doppler  shift  of  3.4  Hz. 


TIME  DELAY  IN  MILLISECONDS 


Figure  84.  Multipath  Intensity  Profile  at  4350  meters  using  DSSS  signal  in  baseband 
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Figure  85.  Multipath  Intensity  Profile  at  4350  meters  -  Bellhop  theoretical  estimate 
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SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  4350  METERS 


Figure  86.  Estimated  Scattering  function  of  the  channel  at  distance  of  4350  meters 


Table  10. 


PATH  WITH  TIME 

DOPPLER 

DOPPLER 

DELAY  (in  msec): 

SHIFT  (HZ) 

SPREAD  (HZ) 

0 

-3.15 

0.66 

1.125 

-3 

1.30 

2.375 

-3 

2.06 

Doppler  spreads  and  shifts  of  the  dominant  paths  for  distance  of  4350  meters 


8.  Received  Signal  at  a  Distance  of  5000  Meters 

The  signal  received  at  time  0113  and  at  a  distance  of  5000  meters  is  illustrated  in 
Figure  87.  This  is  the  most  interesting  case  of  all  since  we  have  a  very  sudden  and  drastic 
reduction  of  the  signal  strength.  The  distance  increased  by  650  meters  and  the  signal- to- 
noise  ratio  reduced  by  9  dB  (a  factor  of  8!).  The  SNR  in  this  case  is  16  dB,  easily  visible 
by  comparing  Figure  87  with  Figure  79.  More  likely,  the  dominant  path  is  beginning  to 
be  unsupported  by  the  channel  geometry.  The  eigenray  propagation  modeled  by  the  Bell¬ 
hop  program  for  the  distance  of  5000  meters  is  presented  in  Figure  88.  The  plot  is  very 
crowded  with  all  the  eigenrays  in  this  distance,  so  we  can  not  make  any  conclusion  for 
the  propagation  with  respect  to  the  drastic  drop  of  the  SNR. 
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Figure  87.  Received  signal  at  distance  of  5000  meters 
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Figure  88.  Eigenrays  plot  for  distance  of  5000  meters 


IN*  use  OF  IFM  chop* 


Figure  89.  Estimated  impulse  response  at  5000  meters 

using  LFM  chirp 
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Figure  90.  Estimated  impulse  response  at  5000  meters  using  DSSS  signal  in  bandpass 
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Figure  91.  Estimated  impulse  response  at  5000  meters  using  DSSS  signal  in  baseband 
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The  multipath  intensity  profile  of  the  channel  at  5000  meters  distance  is  presented 
in  Figure  92.  The  multipath  effect  at  this  distance  is  by  far  the  most  severe  one  of  all  pre¬ 
vious  cases.  We  have  a  total  number  of  six  or  seven  dominant  paths  with  characteristi¬ 
cally  strong  intensities  (in  comparison  with  the  first  path).  The  resulting  multipath  spread 
has  a  value  of  about  6  milliseconds.  Also,  the  dominant  path  in  the  multipath  intensity 
profile  has  a  value  of  7.2  x  10“7 ,  while  in  the  previous  case  it  had  a  value  of  2.6  x  10“5 ,  so 
the  dominant  path’s  intensity  is  about  35  times  weaker  than  previously.  The  Bellhop 
model,  shown  in  Figure  93,  indicates  the  severe  multipath  as  well.  The  match  is  not  per¬ 
fect  but  we  can  conclude  that  the  fading  in  this  distance  is  much  stronger  than  in  4350 
meters  distance.  In  Figure  94  the  scattering  function  of  the  channel  is  displayed.  The  se¬ 
vere  multipath  is  evident.  In  Table  11,  we  summarize  the  estimated  Doppler  spreads  and 
shifts  of  the  dominant  paths.  The  estimated  Doppler  spreads  vary  from  0.84  to  2.36  Hz. 
The  Doppler  shift  of  the  first  path  has  a  negative  value  of  3.1  Hz.  By  comparison,  the 
tone  probes  show  the  same  Doppler  shift  of  -3.1  Hz. 


MULTIPATH  IIMTEMSIT’T  PROFILE  AT  DISTANCE  &□□□  METERS 
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Figure  92.  Multipath  Intensity  Profile  at  5000  meters  using  DSSS  signal  in  baseband 


Figure  93. 


Multipath  Intensity  Profile  at  5000  meters  -  Bellhop  theoretical  estimate 
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SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  5000  METERS 
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Figure  94.  Estimated  Scattering  function  of  the  channel  at  distance  of  5000  meters 


PATH  WITH  TIME 

DOPPLER 

DOPPLER 

DELAY  (in  msec): 

SHIFT  (HZ) 

SPREAD  (HZ) 

0 

-3.1 

0.84 

0.875 

-2.95 

1.51 

1.25 

-2.7 

1.83 

1.75 

-2.8 

1.84 

4.375 

-2.9 

1.7 

5.875 

-2.6 

2.36 

Table  11.  Doppler  spreads  and  shifts  of  the  dominant  paths  for  distance  of  5000  meters 


9.  Received  Signal  at  a  Distance  of  5750  Meters 

The  signal  received  at  time  0143  and  at  a  distance  of  5750  meters  is  illustrated  in 
Figure  95.  In  this  case  the  received  signal  power  reduces  even  more  resulting  in  a  lower 
signal-to-noise  ratio  of  13  dB.  The  plots  of  the  impulse  response  function  are  presented  in 
Figures  96,  97  and  98. 
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Figure  95.  Received  signal  at  distance  of  5750  meters 


Figure  96.  Estimated  impulse  response  at  5750  meters 

using  LFM  chirp 
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Figure  97.  Estimated  impulse  response  at  5750  meters  using  DSSS  signal  in  bandpass 
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|h|2  USE  OF  DSSS  SIGNAL  IN  BASEBAND 


Figure  98.  Estimated  impulse  response  at  5750  meters  using  DSSS  signal  in  baseband 

The  multipath  intensity  profile  at  5750  meters  is  illustrated  in  Figure  99.  In  this 
distance,  we  can  distinguish  six  dominant  paths  in  the  underwater  channel,  whereas  the 
multipath  spread  is  about  8  milliseconds.  Figure  100  shows  the  scattering  function  of  the 
channel.  The  multipath  effect  is  not  as  severe  as  in  the  last  case.  In  Table  12,  the  esti¬ 
mated  Doppler  spreads  and  shifts  are  presented  for  the  dominant  paths.  The  estimated 
Doppler  spread  varies  from  0.94  Hz  to  3.14Hz.  We  estimated  that  the  first  path  has  a 
negative  Doppler  shift  with  a  value  of  3.3  Hz.  By  comparison,  the  tone  probe  shows  a 
negative  Doppler  shift  of  3.5  Hz. 
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Figure  99.  Multipath  Intensity  Profile  at  5750  meters  using  DSSS  signal  in  baseband 
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SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  5750  METERS 
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Figure  100.  Estimated  Scattering  function  of  the  channel  at  distance  of  5750  meters 


PATH  WITH  TIME 

DOPPLER 

DOPPLER 

DELAY  (in  msec): 

SHIFT  (HZ) 

SPREAD  (HZ) 

0 

-3.3 

0.94 

0.75 

-3.15 

1.57 

1.625 

-3.15 

1.50 

3.75 

-2.55 

2.52 

4.625 

-2.95 

1.91 

7.75 

-2.3 

3.14 

Table  12.  Doppler  spreads  and  shifts  of  the  dominant  paths  for  distance  of  5750  meters 

10.  Received  Signal  at  a  Distance  of  6550  Meters 

The  signal  received  at  time  0213  and  at  a  distance  of  6550  meters  is  illustrated  in 
Figure  101.  At  this  distance,  we  are  very  close  to  losing  the  signal  inside  the  noise. 
Probably  that  distance  would  be  the  communication  limit  of  this  underwater  environment. 
The  signal-to-noise  ratio  at  the  distance  of  6.5  kilometers  is  about  7  dB.  The  plots  of  the 
impulse  response  function  for  this  case  are  presented  in  Figures  102,  103  and  104. 
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Figure  101.  Received  signal  at  distance  of  6550  meters 


|h|2  USE  OF  LFM  chirps 


Figure  102.  Estimated  impulse  response  at  6550  meters  using  LFM  chirp 
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Figure  103.  Estimated  impulse  response  at  6550  meters  using  DSSS  signal  in  bandpass 


|h|2  USE  OF  DSSS  SIGNAL  IN  BASEBAND 


Figure  104.  Estimated  impulse  response  at  6550  meters  using  DSSS  signal  in  baseband 

The  multipath  intensity  profde  at  6550  meters  is  illustrated  in  Figure  105.  At  this 
distance,  we  can  distinguish  three  strong  dominant  paths  in  the  underwater  channel,  creat¬ 
ing  a  severe  multipath  environment.  The  multipath  spread  is  small,  in  the  order  of  1.3 
milliseconds.  Figure  106  plots  the  scattering  function  of  the  channel,  revealing  three  large 
components.  In  Table  13  the  estimated  Doppler  spreads  and  shifts  are  presented  for  the 
dominant  paths.  The  estimated  Doppler  spread  varies  from  1 .42  Hz  to  1 .72  Hz.  We  esti¬ 
mated  that  the  first  path  has  a  negative  Doppler  shift  with  a  value  of  4. 1  Hz.  By  compari¬ 
son,  the  tone  probe  experiences  a  negative  Doppler  shift  of  4  Hz. 
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What  is  very  important  to  mention  for  our  research  is  that  even  while  the  SNR  is 
very  low  and  we  may  not  have  a  reliable  communication  channel,  the  information  about 
the  underwater  channel  remains  unambiguous  and  is  almost  unaffected  by  the  relatively 
high  environmental  noise! 


MULTIPATH  INTENSITY  PROFILE  AT  DISTANCE  6550  METERS 


Figure  105.  Multipath  Intensity  Profile  at  6550  meters  using  DSSS  signal  in  baseband 


SCATTERING  FUNCTION  OF  THE  CHANNEL  FOR  DISTANCE  6550  METERS 


Figure  106.  Estimated  Scattering  function  of  the  channel  at  distance  of  6550  meters 
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Table  13. 


PATH  WITH  TIME 

DOPPLER 

DOPPLER 

DELAY  (in  msec): 

SHIFT  (HZ) 

SPREAD  (HZ) 

0 

-4.1 

1.42 

0.625 

-4 

1.73 

1.25 

-4 

1.72 

Doppler  spreads  and  shifts  of  the  dominant  paths  for  distance  of  6550  meters 


11.  Summary  of  the  Results 

Figure  107  presents  a  summary  of  the  multipath  intensity  profiles  for  the  ten  dif¬ 
ferent  cases  of  varying  distance.  Figure  108  presents  a  summary  of  the  channel  scattering 
functions  for  the  ten  different  ranges. 


ESTIMATED  IMPULSE  RESPONSE  FUNCTIONS  FOR  THE  10 
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Figure  107.  Summary  of  the  MIPs  for  the  10  different  cases 
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ESTIMATED  SCATERING  FUNCTIONS  FOR  THE  10  DIFFERENT 

INCIDENT  RANGES 
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Figure  108.  Summary  of  the  scattering  functions  for  the  10  different  cases 


D.  CHAPTER  SUMMARY 

In  this  chapter,  we  processed  the  data  from  an  actual  experiment  and,  using  the 
code  developed  in  the  previous  chapter,  we  determined  the  characteristics  of  the  channel. 
We  described  also  the  conditions  of  this  experiment,  the  format  of  the  data  used  and  also 
referred  in  the  problems  encountered  in  extracting  the  scattering  function  from  the  set  of 
real  data.  The  next  chapter  presents  a  summary  with  conclusions  and  goals  achieved  in 
this  research. 
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VI.  CONCLUSIONS  AND  FUTURE  WORK 


A.  CONCLUSIONS 

This  thesis  examined  the  behavior  of  the  underwater  channel  in  which  the  Seaweb 
communication  network  operates.  After  developing  the  theory  describing  the  underwater 
channel  from  both  physics  and  electrical  engineering  perspectives,  and  discussing  how 
the  channel  characteristics  affect  the  acoustic  communication  signal,  we  described  an  ef¬ 
ficient  method  of  measuring  those  characteristics,  even  in  the  presence  of  noise. 

The  next  step  was  to  simulate  a  time-varying  underwater  channel  with  multipath 
characteristics.  After  passing  a  bandpass  communication  signal  through  this  channel,  we 
applied  the  method  developed  in  the  previous  chapter  to  measure  the  characteristics  of 
the  artificial  channel.  The  purpose  of  this  action  was  twofold:  first,  the  simulated  channel 
served  as  a  test  channel  to  ensure  that  the  method  works  as  expected  and  secondly  to  ob¬ 
serve  the  influence  of  the  gradually  increased  noise  levels  on  the  results  of  the  estimation. 

The  results  were  quite  satisfactory  since  they  were  accurate  even  for  signal  to 
noise  ratio  of  7  dB.  Finally,  we  described  an  actual  ocean  experiment  and  the  various 
types  of  signals  sent  as  channel  probes.  Each  signal  type  served  a  different  purpose  with 
respect  to  determining  the  effect  of  the  channel  on  the  acoustic  communication  signal. 
The  different  types  of  signal  sent  gave  us  the  opportunity  to  examine  other  ways  of  chan¬ 
nel  estimation  besides  the  one  we  described  in  Chapter  IV. 

The  key  findings  from  this  work  are  as  follows: 

•  LFM  chirps  and  the  matched  filter  method  gave  us  a  very  clear  picture  for 
the  estimated  impulse  response  of  the  channel.  However,  the  low  meas¬ 
urement  rate  of  this  method  led  to  great  aliasing  in  the  scattering  function, 
which  made  Doppler  estimation  impossible. 

•  The  DSSS  method  worked  very  well  in  both  the  bandpass  and  baseband 
cases,  so  that  the  transition  from  bandpass  to  baseband  prior  to  channel  es¬ 
timation  is  unnecessary. 
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•  The  results  of  the  DSSS  method  gave  a  noisy  picture  for  the  impulse  re¬ 
sponse  because  of  the  format  of  the  DSSS  signal.  The  multipath  intensity 
profile  on  the  other  hand  gave  a  good  estimate  of  the  channel  (still  noisy). 
The  estimate  of  the  scattering  function  from  the  DSSS  method  came  out 
very  clear  and  very  accurate,  giving  us  unambiguous  information  for  the 
channel. 

•  Each  of  the  different  paths  in  the  multipath  has  a  different  value  of  Dop¬ 
pler  shift  and  Doppler  spread.  The  weaker  components  of  the  impulse  re¬ 
sponse  usually  have  higher  Doppler  spreads  than  the  stronger  ones,  and 
the  impulse  response  component  corresponding  to  the  first  path  usually 
has  the  highest  Doppler  shift. 

•  The  tonal  comb  signal  gave  us  a  value  of  Doppler  shift  very  close  or  ex¬ 
actly  the  same  as  that  of  the  impulse  response  component  corresponding  to 
the  first  path. 


B.  FUTURE  WORK 

Since  we  have  shown  the  feasibility  of  measuring  the  important  characteristics  of 
the  underwater  channel,  the  next  step  is  to  implement  the  DSSS  signal  and  the  appropri¬ 
ate  algorithm  inside  the  Seaweb  modem.  A  natural  implementation  is  to  send  a  known 
DSSS  signal  in  the  beginning  of  the  communication  (inside  the  RTS  transmission). 

A  separate  important  task  will  be  associating  the  underwater  channel  types  with 
the  available  communication  schemes  and  their  parameters  such  as  frequency  band  of  op¬ 
eration,  modulation  scheme  (coherent,  non  coherent  and  specific  type),  data  rate,  modem 
output  power,  and  error  correction  coding  type,  which  will  give  us  the  minimum  prob¬ 
ability  of  error.  The  match  up  can  be  done  by  researching  the  several  experiments  that 
have  been  done  in  different  ocean  environments,  by  organizing  new  ones  towards  this 
specific  goal,  and  by  simulating  the  communication  system.  We  anticipate  that  there  must 
be  a  limiting  case  of  severity  in  the  channel  conditions  beyond  which  we  cannot  use  any 
phase  coherent  techniques  with  equalization  and  we  are  compelled  to  use  the  more  robust 
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non-coherent  FSK,  since  this  appears  to  be  less  susceptible  to  channel  impairments  and 
seems  to  provide  robust  communications  and  relatively  high  data  rates. 

The  ultimate  scope  of  future  work  is  to  build  an  adaptive  modem  that  will  esti¬ 
mate  the  character  of  the  channel  using  the  DSSS  method  and,  based  on  that,  will  decide 
on  the  communication  parameters  which  optimize  communication  between  neighboring 
pairs  of  Seaweb  nodes. 
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APPENDIX. 


MATLAB  CODES 


In  this  appendix,  we  present  the  Matlab  codes  used  in  this  work  to  generate  the  re¬ 
sults  and  the  plots  for  the  artificial  and  the  actual  channel. 


%  GENERATION  OF  DSSS  (GOLD  CODE  WITH  LENGTH  2047) 
%  AND  PROBE  SIGNAL  WITH  LENGTH  8  bits 
clc 
clear 
close  all 
format  long 

load  gold2047_l 
codel=code(l,:); 
code2=code(2,:); 
s  1  =sign(randn(  1,8)); 
sa=(codel'*sl)+(i*code2'*sl); 
s=reshape(sa,  1 ,2047 *length(s  1 )) ; 
s=s+0. 1  *randn(  1 ,2047*length(s  1 )); 


%  UPSAMPLING  THE  DATA 

%  sampleperchip  =  ones(l,l); 

%  hip  =  sampleperchip'*s2; 

%  s  =  reshape(hlp,  1  ,prod(size(hlp))); 
s=s'; 


%  CREATE  DATA  FOR  USE 

si=D; 

for  i=  1:201 
ss(i,i)=0; 
sl=[sl  s(i,  1)]; 
ss(i,:)=sl; 

end 

ss=conj(ss); 

sa=[]; 

for  i=l:(200) 

sa=[sa ;  s(l:200)']; 

end 
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for  i=l  :(length(s)-200) 
i 

sb(i,:)=s(i:i+200)'; 

end 


%  GENERATION  OF  IMPULSE  RESPONSE  CORRESPONDING 
%  TO  THE  DIRECT  PATH  OF  THE  ARTIFICIAL  CHANNEL 


clc 
clear 
close  all 
format  long 

k=l  :(8*2047+200); 

hl=(0.1*sin(0.0004*k-0.2)-0.16*sin(0.00032*k+0.7)+0.05*cos(0.0006*k+0.3)- 
0.13*sin(0.00049*k-0.2)+0.07*sin(0.00064*k-0.5)+0.11*sin(0.00057*k+0.1)- 
0.07*cos(0.00055*k)+0.16*sin(0.00052*k+0.6)-0.06*cos(0.00048*k- 
0.1)+0.05*cos(0.00045*k+0.9)-0.03*sin(0.00041*k-0.4))*0.6; 
h  1  =h  1  -mean(h  1 )+ 1 ; 

%  UPSAMPLING  THE  IMPULSE  RESPONSE 

sampleperchip  =  ones(l,l); 
hip  =  sampleperchip'*hl; 
h  =  reshape(hlp,  1  ,prod(size(hlp))); 
h=h'; 


%  PLOTTING  IMPULSE  RESPONSE  MAGNITUDE  AND  PHASE 
%  AS  A  FUNCTION  OF  ABSOLUTE  TIME 


t=0:0.25:((8*2047-l+200)*0.25); 

figure(l) 

plot(t,abs(h)) 

xlabel('TIME  IN  msec') 

ylabel('MAGNITUDE  OF  IMPULSE  RESPONSE') 

axis([0  max(t)  0  2]) 

grid 

figure(2) 

plot(t,phase(h)) 

xlabel('TIME  IN  msec') 

ylabel('PHASE  OF  IMPULSE  RESPONSE') 

grid 


96 


figure(3) 

H=fft(h-mean(h),4096); 

H=fftshift(H); 

Phh  =  H.*conj(H)  /  4096; 
f  =  4000*(0:4095)/4096; 
plot(f,(Phh)) 


%  PLOTTING  THE  COHERENCE  FUNCTION  OF 
%  THIS  IMPULSE  RESPONSE 

rho=xcorr(h,'coeff); 

dt=-0.25*33 150/2:0.25:0.25*33 150/2; 

figure(4) 

plot(  1 0A-3  *dt,(rho)) 

xlabel('  dt  in  seconds  ') 

ylabel('  Normalised  autocorrelation ') 

title('  COHERENCE  FUNCTION  OF  THE  DIRECT  PATH  ') 
grid 
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%  GENERATION  OF  IMPULSE  RESPONSE  CORRE 
%  SPONDING  TO  THE  SECOND 
%  PATH  OF  THE  ARTIFICIAL  CHANNEL 


clc 
clear 
close  all 
format  long 


k=l:(2047*8+200); 

hl=(-0.2*sin(0.00035*k-0.2)-0.1*sin(0.00036*k+0.9)+0.09*cos(0.000411*k+0.6)- 
0.09*sin(0.000361*k-0.4)+0.13*sin(0.00047*k-.5)+0.15*sin(0.00052*k+0.3)- 
0.13*cos(0.00048*k-0.6)-0.07*sin(0.00057*k+0.7)-0.06*cos(0.00046*k)- 
0.04*cos(0.00043*k-0.6)-0.04*sin(0.00037*k-0.4))*1.3; 
h  1  =h  1  -mean(h  1  )+0 . 5 ; 

%  UPSAMPLING  THE  IMPULSE  RESPONSE 

sampleperchip  =  ones(l,l); 
hip  =  sampleperchip'*hl; 
h  =  reshape(hlp,  1  ,prod(size(hlp))); 
ha=h'; 


%  PLOTTING  IMPULSE  RESPONSE  MAGNITUDE  AND 
%  PHASE  AS  A  FUNCTION  OF  ABSOLUTE  TIME 


t=0:0.25:((2047*8-l+200)*0.25); 

figure(l) 

plot(t,abs(ha)) 

xlabel('TIME  IN  msec’) 

ylabel('MAGNITUDE  OF  IMPULSE  RESPONSE') 

axis([0  max(t)  0  2]) 

grid 

figure(2) 

plot(t,phase(ha)) 

xlabel('TIME  IN  msec') 

ylabel('PHASE  OF  IMPULSE  RESPONSE') 

grid 

figure(3) 

H=fft(ha-mean(ha),4096); 

H=fftshift(H); 

Phh  =  H.*conj(H)  /  4096; 


98 


f  =  4000*(0:4095)/4096; 
plot(f,(Phh)) 


%  PLOTTING  THE  COHERENCE  FUNCTION  OF  THIS 
%  IMPULSE  RESPONSE 

rho=xcorr(ha,'coeff); 

dt=-0.25*33 150/2:0.25:0.25*33 150/2; 

figure(4) 

plot(  1 0A-3  *dt,(rho)) 

xlabel('  dt  in  seconds  ') 

ylabel('  Normalized  autocorrelation ') 

title('  COHERENCE  FUNCTION  OF  THE  SECOND  PATH  (TIME  DELAY  10  msecs)  ’) 
grid 
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%  GENERATION  OF  IMPULSE  RESPONSE  CORRE- 
%  SPONDING  TO  THE  THIRD  PATH  OF  THE  ARTIFICIAL 
%  CHANNEL 


clc 
clear 
close  all 
format  long 

k=l:(2047*8+200); 

hl=(0.1*sin(0.00030*k-0.4)- 

0.4*sin(0.00046*k+0.9)+0.07*cos(0.00028*k+0.3)+0.07*sin(0.00041*k-0.1)- 
0.14*sin(0.00035*k-.5)-0.18*sin(0.00026*k+0.1)+0.14*cos(0.00035*k-0.6)- 
0.08*sin(0.00037*k+0.7)+0.07*cos(0.00034*k-l.l)-0.09*cos(0.00029*k-0.5)- 
0.08*sin(0.00036*k-0.2))*0.5; 
h  1  =h  1  -mean(h  1  )+0 .2 ; 

%  UPSAMPLING  THE  IMPULSE  RESPONSE 

sampleperchip  =  ones(l,l); 
hip  =  sampleperchip'*hl; 
h  =  reshape(hlp,  1  ,prod(size(hlp))); 
hb  IV; 


%  PLOTTING  IMPULSE  RESPONSE  MAGNITUDE  AND 
%  PHASE  AS  A  FUNCTION  OF  ABSOLUTE  TIME 

t=0:0.25:((2047*8-l+200)*0.25); 

figure(l) 

plot(t,abs(hb)) 

xlabel('TIME  IN  msec’) 

ylabel('MAGNITUDE  OF  IMPULSE  RESPONSE') 
title('IMPULSE  RESPONSE  AT  TIME  DELAY  10  MILLISECONDS') 
axis([0  max(t)  0  1]) 
grid 

figure(2) 

plot(t,phase(hb)) 

xlabel('TIME  IN  msec') 

ylabel('PHASE  OF  IMPULSE  RESPONSE') 

grid 

figure(3) 

H=fft(hb-mean(hb)  ,4096) ; 

H=fftshift(H); 
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Phh  =  H.*conj(H)  /  4096; 
f  =  4000*(0:4095)/4096; 
plot(f,(Phh)) 


%  PLOTTING  THE  COHERENCE  FUNCTION  OF  THIS 
%  IMPULSE  RESPONSE 


rho=xcorr(hb,'coeff); 

dt=-0.25*33 150/2:0.25:0.25*33 150/2; 

figure(4) 

plot(  1 0A-3  *dt,(rho)) 

xlabel('  dt  in  seconds  ') 

ylabel('  Normalised  autocorrelation ') 

title('  COHERENCE  FUNCTION  CORRESPONDING  TO  THE  THIRD  PATH  (TIME 

DELAY  50  msecs) ') 

grid 


101 


%  COMPARISON  OF  THE  THREE  COMPONENTS  OF  THE  IM 
%  PULSE  RESPONSE  OF  THE  ARTIFICIAL  CHANNEL 


clc 
clear 
close  all 
format  long 

k=l  :(2047*8+200); 


hl=(0.1*sin(0.00030*k-0.4)- 

0.4*sin(0.00046*k+0.9)+0.07*cos(0.00028*k+0.3)+0.07*sin(0.00041*k-0.1)- 

0.14*sin(0.00035*k-.5)-0.18*sin(0.00026*k+0.1)+0.14*cos(0.00035*k-0.6)- 

0.08*sin(0.00037*k+0.7)+0.07*cos(0.00034*k-l.l)-0.09*cos(0.00029*k-0.5)- 

0.08*sin(0.00036*k-0.2))*0.5; 

hb=h  1  -mean(h  1  )+0 .2 ; 


hl=(-0.2*sin(0.00035*k-0.2)-0.1*sin(0.00036*k+0.9)+0.09*cos(0.000411*k+0.6)- 
0.09*sin(0. 000361  *k-0.4)+0.13*sin(0.00047*k-.5)+0.15*sin(0.00052*k+0. 3)- 
0.13*cos(0.00048*k-0.6)-0.07*sin(0.00057*k+0.7)-0.06*cos(0.00046*k)- 
0.04*cos(0.00043*k-0.6)-0.04*sin(0.00037*k-0.4))*  1 .3; 
ha=h  1  -mean(h  1  )+0 . 5 ; 


hl=(0.1*sin(0.0004*k-0.2)-0.16*sin(0.00032*k+0.7)+0.05*cos(0.0006*k+0.3)- 
0. 13*sin(0.00049*k-0.2)+0.07*sin(0.00064!,:k-0.5)+0. 1 1  *sin(0.00057*k+0. 1)- 
0.07*cos(0.00055*k)+0.16*sin(0.00052*k+0.6)-0.06*cos(0.00048*k- 
0.1)+0.05*cos(0.00045*k+0.9)-0.03*sin(0.00041*k-0.4))*0.6; 
h=h  1  -mean(h  1 )+ 1 ; 


t=0:0.00025:((2047*8-l+200)*0.25/1000); 

figure(l) 

subplot(3,l,l) 

plot(t,abs(h)) 

title('IMPULSE  RESPONSE  AT  TIME  DELAY  ZERO') 

axis([0  max(t)  0  2]) 

grid 


subplot(3,l,2) 

plot(t,abs(ha)) 

title('IMPULSE  RESPONSE  AT  TIME  DELAY  10  MILLISECONDS') 
ylabel('MAGNITUDE  OF  IMPULSE  RESPONSE') 
axis([0  max(t)  0  2]) 
grid 

subplot(3,l,3) 
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plot(t,abs(hb)) 

xlabel('TIME  IN  SECONDS') 

title('IMPULSE  RESPONSE  AT  TIME  DELAY  50  MILLISECONDS') 

axis([0  max(t)  0  1]) 

grid 


htel=[h;  zeros(39, 16576);  ha;  zeros(159, 16576);  hb;  zeros(l,16576)]; 

%  MULTIPATH  INTENSITY  PROFILE  OF  THE  ARTIFICIAL 
%  CHANNEL 


figure(2) 

a=conj  (htel).  *htel; 
phih=0 . 5  *  sum(a')/ 16576; 
plot(0:0. 25:201  *0.25, phih) 
axis([-0. 11  50.5  0  0.53]) 
grid 

title('  MULTIPATH  INTENSITY  PROFILE  (MIP) ') 
xlabel('  TIME  DELAY  IN  MILLISECONDS  ') 
ylabel('  AVERAGE  MULTIPATH  INTENSITY  ') 
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%  TEST  CHANNEL  ESTIMATION 


clc 
clear 
close  all 


load  dataprosepejergasia; 
load  impulse  1; 
load  impulse05; 
load  impulse02; 


%  DATA  GENERATION 

%  FIRST  COMPONENT  OF  IMPULSE  RESPONSE 
%  SECOND  COMPONENT  OF  IMPULSE  RESPONSE 
%  THIRD  COMPONENT  OF  IMPULSE  RESPONSE 


htel=[h.';  zeros(39, 16576);  ha.';  zeros(  159, 16576);  hb.';  ]; 

%  DATA  PASSING  THROUGH  THE  CHANNEL 
%  RECEIVED  SIGNAL  GENERATION 
%  PLUS  AWGN 


s=fliplr(sb); 
ss=s(l ,:); 

spro=zeros(200,201); 
for  i=  1:200 

spro(20 1  -i,  :)=[ss(i+ 1 :20 1 )  zeros(  1  ,i)] ; 
end 

st=s(16176,:); 

smeta=zeros(200,201); 

for  i=  1:200 

smeta(i,:)=[zeros(  1  ,i)  st(  1:201  -i)] ; 
end 

stel=[  spro  ;  s;  smeta]; 

for  i=l:  16576 

rcv(i)=stel(i, : )  *htel( :  ,i) ; 
end 

rcv=rcv+0 .44  *  (randn(  1 , 1 6 5 76)+j  *randn(  1,16576)); 
figure(l) 

plot(real(rcv),imag(rcv),'*') 


for  i=l:8 

c(i)=sb(2047 *  (i- 1 )+ 1 ) ; 
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end 

sl=-c./(l+j); 

load  gold2047_l 

codel=code(l,:); 

code2=code(2,:); 

sa=(code  l'*s  1  )+(j  *code2'*s  1 ); 

s2=reshape(sa,  1 ,2047 *length(s  1 )) ; 

%  UPSAMPLING  THE  DATA 


sampleperchip  =  ones(l,l); 

hip  =  sampleperchip'*s2; 

sp  =  reshape(hlp,  1  ,prod(size(hlp))); 

sp=sp.'; 


%  CHANNEL  ESTIMATION 


sc=conj(sp); 

d=0; 

for  g=210: 100: 14300 
g 

d=d+ 1 ; 
for  k=  1:205 
sum  1=0; 
for  m=  1:2047 

boh=rc  v(g+m)  *  sc(g+m-k+ 1 ) ; 
sum  1  =boh+sum  1 ; 
end 

hestim(d,k)=sum 1/(2047*2); 
end 
end 


%  PLOTTING  THE  ESTIMATE  OF  IMPULSE 
%  RESPONSE  AS  WELL  AS  THE  ACTUAL  IMPULSE 
%  RESPONSE  FOR  COMPARISON  PURPOSES 


t=0:0.0025:0.0025*  16575; 

tl=210*0.00025:0. 025: 14300*0.00025; 

figure(2) 

subplot(2,l,l) 

plot(t,h) 

title('  ACTUAL  IMPULSE  RESPONSE  OF  THE  DIRECT  PATH  ') 
ylabel(’  MAGNITUDE  OF  IMPULSE  RESPONSE  ') 
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grid 

hesta=abs(hestim(:,  1)); 

subplot(2, 1 ,2) 
plot(tl,hesta) 

title('  ESTIMATED  IMPULSE  RESPONSE  OF  THE  DIRECT  PATH  ') 

xlabel('  ABSOLUTE  TIME  IN  SECONDS  ’) 

ylabel(’  MAGNITUDE  OF  IMPULSE  RESPONSE  ’) 

axis([0  4.5  0.8  1.2]) 

grid 

figure(3) 

subplot(2,l,l) 

plot(t,ha) 

title(’  ACTUAL  IMPULSE  RESPONSE  OF  THE  SECOND  PATH  (TIME  DELAY  10 
mecs) ') 

ylabel('  MAGNITUDE  OF  IMPULSE  RESPONSE  ’) 
grid 

hestb=hestim( :  ,4 1 ) ; 

subplot(2, 1 ,2) 
plot(tl,hestb) 

title('  ACTUAL  IMPULSE  RESPONSE  OF  THE  SECOND  PATH  (TIME  DELAY  10 
mecs) ') 

xlabel(’  ABSOLUTE  TIME  IN  SECONDS  ') 
ylabel('  MAGNITUDE  OF  IMPULSE  RESPONSE  ') 
axis([0  4.5  0  1]) 
grid 

figure(4) 

subplot(2,l,l) 

plot(t,hb) 

title('  ACTUAL  IMPULSE  RESPONSE  OF  THE  THIRD  PATH  (TIME  DELAY  50 
mecs) ') 

ylabel(’  MAGNITUDE  OF  IMPULSE  RESPONSE  ') 
grid 

hestc=abs(hestim(:  ,20 1 )); 

subplot(2,l,2) 

plot(tl,hestc) 

title('  ACTUAL  IMPULSE  RESPONSE  OF  THE  THIRD  PATH  (TIME  DELAY  50 
mecs) ') 
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xlabelC  ABSOLUTE  TIME  IN  SECONDS  ') 
ylabel('  MAGNITUDE  OF  IMPULSE  RESPONSE  ') 
axis([0  4.5  0.05  0.4]) 
grid 


%  PLOTTING  ESTIMATE  OF  MIP 


a=conj  (hestim).  *hestim; 
phih=0 . 5  *  sum(a)/ 1380; 
figure(8) 

plot(0:0. 25:204*0.25, phih) 
axis([-0.1  51  0  0.55]) 

title('  ESTIMATED  MULTIPATH  INTENSITY  PROFILE  OF  ARTIFICIAL  UNDER¬ 
WATER  CHANNEL ') 

xlabelC  TIME  DELAY  IN  MILLISECONDS  ') 
ylabel('  AVERAGE  MULTIPATH  INTENSITY  ') 
grid 


%  ESTIMATING  SCATTERING  FUNCTION 
%  OF  ARTIFICIAL  CHANNEL 


figure(9) 
for  k=  1:205 

phihtau(:,k)=xcorr(hestim(:,k)-mean(hestim(:,k))); 

end 

for  k=  1:205 

phihtau  1  ( :  ,k)=phihtau( :  ,k) ; 
end 

for  k=  1:205 

Sfin( :  ,k)=fft(phihtau  1  ( :  ,k)) ; 
Sfm(:,k)=fftshift(Sfm(:,k)); 
end 

f  =  40*(-140:140)/281; 
fl=f; 


%  ESTIMATING  DOPPLER  SPREAD 

gn=(abs(Sfin(:,l)))./sum(abs(Sfin(:,l))); 

Bdl=sqrt(sum(fl  ,A2.*gn)-(sum(fl  ,*gn)).A2); 

gn=(abs(Sfin(:  ,4 1 )))./ sum(abs(Sfin( :  ,4 1 ))) ; 

Bd2=sqrt(sum(fl  ,A2.*gn)-(sum(fl  ,*gn)).A2); 
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gn=(abs(Sfin(:,201)))./sum(abs(Sfin(:,201))); 

Bd3=sqrt(sum(fl  ,A2.*gn)-(sum(fl  ,*gn)).A2); 

Bdl2=sqrt(sum(fl.A2.*abs(Sfin(:,l)))./sum(abs(Sfin(:,l)))); 
Bd22=sqrt(sum(fl .  A2 .  *abs(Sfin( :  ,4 1  )))./sum(abs(Sfin( :  ,4 1 )))); 
Bd32=sqrt(sum(fl.A2.*abs(Sfin(:,201)))./sum(abs(Sfin(:,201)))); 

%  PLOTTING  SCATTERING  FUNCTION 
%  OF  ARTIFICIAL  CHANNEL 


tau=0:0.25:204*0.25; 

[X,Y]  =  meshgrid(tau,f); 

Sfinmag=abs(Sfin); 

Sfinmagnorm=Sfinmag./max(max(Sfinmag)); 

plot3(X,Y,Sfinmagnorm) 

title('  ESTIMATED  (NORMALIZED)  SCATTERING  FUNCTION  OF  THE  CHAN¬ 
NEL  ’) 

xlabelC  TIME  DELAY  IN  mseconds  ’) 
ylabelC  FREQUENCY  IN  Hz  ') 
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%  IMPULSE  RESPONSE  PROFILE  ESTIMATION  USING  LFM  CHIRP 


clc 
clear 
close  all 
format  short 

SENT=wavread('Probe_PCM'); 
RCVD=wavread('D2 1 0 1 43  AN’); 
SENT=SENT(48000:50400,1); 
RCVD=RCVD(850354: 1450000); 

d=0; 

for  g=l  :6:(length(RCYD)-2402) 
g 

d=d+l; 
sum  1=0; 

for  m=  1:2401 

boh=SENT(m)*(RCVD(m+g-l)); 
sum  1  =boh+sum  1 ; 
end 

hestim(d)=sum 1/(2401); 
end 

A=abs(hestim); 

%  for  n=l  :length(hestim) 

%  if  abs(A(n))<0.08*max(A) 

%  A(n)=0; 

%  end 
%  end 

for  n=l:40 

hestimate(n, :  )=hestim((n- 1  )*2000+ 1 :2000*n); 
end 

figure(l) 

time=0:0.25:39*0.25; 
tau=0:0.25/2000:399*0.25/2000; 
[X,Y]=meshgrid(tau,time); 
waterfall(X,Y,(abs(hestimate(:,  1:400))).  A2) 
title('|h|A2  USE  OF  LFM  chirps  ') 
xlabel(  time  delay ') 
ylabel(  absolute  time ') 
grid 

colorbar 
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%  IMPULSE  RESPONSE  PROFILE  ESTIMATION  USING  DSSS  SIGNAL  IN 
BANDPASS 


clc 
clear 
close  all 
format  short 

SENT=wavread('junk_2'); 

RCVD=wavread('D2 1011 3  AN'); 

RCVD=RCVD(3630001 :5630000); 

R=fft(RCVD); 

R=fftshift(R); 

a=-3630001+5630000; 

R(a/2- 1  :a/2+2)=0; 

R=ifftshift(R); 

RCVD=ifft(R); 
clear  R 

d=0; 

for  g=75000:4094*3: 1900000 
g 

d=d+l; 

for  k=l:6: 1501 
sum  1=0; 
for  m=  1:4094 

boh=RCVD(g+m)*(SENT(g+31420-68250+m-k+l)); 
sum  1  =boh+sum  1 ; 
end 

hestim(d,(k+5)/6)=suml/(4094); 

end 

end 

a=size(hestim); 

figure(9) 

time=0:2047*3*0.5/12000:(a(l)-l)*2047*3*0. 5/12000; 
tau=0:6/48000:(a(2)-l)*6/48000; 

[X,Y]=meshgrid(tau,time); 
waterfall(X,Y,(abs(hestim)).A2) 
xlabel('  time  delay ') 
ylabel('  absolute  time ') 

title('  |h|A2  USE  OF  DSSS  SIGNAL  IN  BANDPASS  ') 

axis  tight 

colorbar 
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%  IMPULSE  RESPONSE  PROFILE  ESTIMATION  USING  DSSS  SIGNAL  IN  BASE¬ 
BAND 

clc 
clear 
close  all 
format  short 

SENT=wavread('junk_2'); 

RCVD=wavread('D2 1011 3  AN'); 

RCVD=RCVD(3630001 :5630000); 

R=fft(RCVD); 

R=fftshift(R); 
a=-363000 1+5630000; 

R(a/2- 1  :a/2+2)=0; 

R=ifftshift(R); 

RCVD=ifft(R); 
clear  R 

theta=1.03*pi; 

nl=[l:length(SENT)]; 

pl=SENT.*cos(2*pi*12000*nl'/48000); 

Pl=fft(pl); 
clear  pi 

Pl=fftshift(Pl); 

Pl=[zeros(750000,l);  Pl(750001: 1250000,1);  zeros(length(Pl)- 125 0000,1)]; 

Bl=ifftshift(Pl); 

clear  PI 

bl=ifft(Bl); 

clear  B 1 

pl=SENT.*cos(2*pi*12000*nl'/48000); 

Plq=fft(plq); 
clear  plq 
P 1  q=fftshift(P  1  q) ; 

Plq=[zeros(750000,l);Plq(750001: 1250000,1);  zeros(length(Plq)-1250000,l)]; 

B 1  q=ifftshift(P  1  q); 
clear  Plq 
blq=ifft(Blq); 
clear  Blq 

n2=[l  :length(RCVD)]; 

p2=RCVD.*cos(2*pi*12000*n2'/48000+theta); 

P2=fft(p2); 
clear  p2 

P2=fftshift(P2); 
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P2=[zeros(600000, 1 );  P2(60000 1 : 1 400000, 1 );  zeros(length(P2)- 1 400000, 1 )] ; 

B2=ifftshift(P2); 

clear  P2 

b2=ifft(B2); 

clear  B2 

p2q=RCVD.*cos(2*pi*12000*n2'/48000+theta+pi/2); 

P2q=fft(p2q); 
clear  p2q 
P2q=fftshifl(P2q); 

P2q=[zeros(600000, 1 );  P2q(60000 1 : 1400000, 1 );  zeros(length(P2q)- 1 400000, 1 )] ; 

B2q=ifftshift(P2q); 

clear  P2q 

b2q=ifft(B2q); 

clear  B2q 

clear  SENT 

clear  RCVD 

sent=real(b  1  )+j  *real(b  1  q); 
rcvd=real(b2)+j  *real(b2q); 
clear  hestim 

d=0; 

for  g=75000:2047: 1900000 
g 

d=d+l; 

for  k=l:6: 1501 
sum  1=0; 
for  m=  1:4094 

boh=rcvd(g+m)*conj(sent(g+31420-68250+m-k+l)); 
suml=boh+suml ; 
end 

hestim(d,(k+5)/ 6)=sum 1 1(2 *4094) ; 
end 
end 

a=size(hestim); 

figure(l) 

time=0:2047*0.5/(6*4000):(a(l)-l)*2047*0. 5/(6*4000); 
tau=0:6/48000:(a(2)-l)*6/48000; 

[X,Y]=meshgrid(tau,time); 
waterfall(X,Y,(abs(hestim)).A2) 
xlabel('  time  delay ') 
ylabel('  absolute  time ') 
grid 

colorbar 
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%  MULTIPATH  INTENSITY  PROFILE  -  SCATTERING  FUNCTION  - 
%  DOPPLER  SPREADS  AND  SHIFTS  OF  THE  PATHS  (DSSS  BASEBAND) 

clc 
clear 
close  all 
format  short 

load('IMP01 13'); 

a=size(hestim); 

figure(l) 

time=0:2047*0.25/12000:(a(l)-l)*2047*0.25/12000; 

tau=0:6/48000:(a(2)-l)*6/48000; 

[X,Y]=meshgrid(tau,time); 
waterfall(X,Y,(abs(hestim)).A2) 
xlabel('  time  delay ') 
ylabel('  absolute  time ') 

title('  |h|A2  USE  OF  DSSS  SIGNAL  IN  BASEBAND  ') 

axis  tight 

colorbar 

H=hestim.  *coni  (hestim); 
fori=l:150 

[a(i),b(i)]=find(H==max(H(i,  1 : 1 6))); 
end 

for  i=15 1:286 

[a(i),b(i)]=find(H==max(H(i,  1 6 : 30))) ; 
end 

for  i=287:380 

[a(i),b(i)]=find(H==max(H(i, 28:38))); 
end 

for  i=38 1:555 

[a(i),b(i)]=find(H==max(H(i,3 8 : 54))) ; 
end 

for  i=556:800 

[a(i),b(i)]=find(H==max(H(i, 53:75))); 
end 

for  i=80 1:892 

[a(i),b(i)]=find(H==max(H(i,73 : 88))); 
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end 


order=[a;  b]; 
for  i=l:892 

HS(a(i), :)=[H(a(i),(b(i):25 1 ))  zeros(  1  ,b(i))] ; 
end 

HS=HS(:, 1:251); 
figure(2) 

waterfall(X,Y,HS) 
for  i=l:892 

hestimate(a(i),:)=[hestim(a(i),(b(i):251))  zeros(l,b(i))]; 
end 

hestimate=hestimate(:,  1 : 1 69); 
figure(3) 

Hcrit=hestimate.*conj(hestimate); 

MIP=sum(Hcrit)/length(Hcrit) ; 

plot(0 : 0 . 1 2 5 : 0 . 1 25  *  (length(MIP)- 1  ),MIP) 

title('  MULTIPATH  INTENSITY  PROFILE  AT  DISTANCE  5000  METERS  ') 
xlabel('  TIME  DELAY  IN  MILLISECONDS  ') 
ylabel('  INTENSITY ') 
grid 

for  k=l :  169 

phihtau( :  ,k)=xcorr(hestimate( :  ,k)-mean(hestimate( :  ,k))); 
end 

for  k=l :  169 

phihtau  1  ( :  ,k)=phihtau( :  ,k) ; 
end 

for  k=l :  169 

Sfin( :  ,k)=fft(phihtau  1  ( :  ,k)) ; 

Sfm(:,k)=fftshift(Sfm(:,k)); 

end 

q=length(HS); 
for  i=  1 : 1 69 

if  max(abs(Sfin(:,i)))<0.01*max(max(abs(Sfin))) 

Sfin(:,i)=0; 

end 

end 
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figure(4) 
binf=-89 1:891; 

frequency=23. 449/2 *binf/891 ; 
bint=l :  169; 
tau=0.125*bint; 

[X,Y]=meshgrid(tau, frequency); 
waterfall(X’,Y',(abs(Sfm))') 

title('  SCATTERING  FUNCTION  OF  THE  CHANNEE  FOR  DISTANCE  5000  ME¬ 
TERS  ’) 

xlabel('  TIME  DELAY  IN  MILLISECONDS  ') 
ylabel('  FREQUENCY  IN  HZ  ') 
axis  tight 


%  FIND  THE  STRONGEST  PATHS 


for  d=  1:15 

[e,r(d)]=find(abs(Sfin)==max(max(abs(Sfin)))); 

Sband( :  ,d)=Sfin( :  ,r(d)); 

fo(d)=(sum((frequency').*abs(Sfm(:,r(d))))./sum(abs(Sfin(:,r(d))))); 

Bd(d)=sqrt(sum((frequency'-fo(d)).A2.*abs(Sfin(:,r(d))))./sum(abs(Sfin(:,r(d))))); 

Sfm(:,r(d))=0; 

end 


SET=[r'  Bd'  fo']; 
[s,c]=sort(r); 
set=[s',Bd(c)'  fo(c)']; 
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